MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  iUNCAU  OF  STANDARDS- 1»«* “ A 


&D  A1 39656 


ARL-STRUC-REPORT-393 


AR -002-896 


DEPARTMENT  OF  DEFENCE  SUPPORT 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  LABORATORIES 

MELBOURNE.  VICTORIA 

STRUCTURES  REPORT  383 


ON  THE  GENE8I8  OF  RELIABILITY  MODELS 


6.  D.  MALLINSON 


THE  UNITED  STATES  NATIONAL 
TECHNICAL  INFORMATION  8ERVK* 
IS  AUTHORISED  TO  ___ 

REPRODUCE  AND  SELL  THIS  REPORT 


Approvd  for  Public  MImm 


(C) 


MBBMUA  1982 


file  copy  ** 


OS  TB 


:  ©<8. 


JULY  1982 


AR-OQ2-806 


DEPARTMENT  OF  DEFENCE  SUPPORT 
DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  LABORATORIES 


STRUCTURES  REPORT  393 


ON  THE  GENESIS  OF  RELIABILITY  MODELS 


by 

G.  O.  MALUNSON 


SUMMARY 

A  method  for  the  derivation  of  reliability  models  for  fatigue  is  presented  and  its 
application  demonstrated. 


DTIC 


©  COMMONWEALTH  OF  AUSTRALIA  1982 

POSTAL  ADDRESS:  Director,  Aeronautical  Research  Laboratories, 

Box  4331,  P.O.,  Melbourne,  Victoria,  3001,  Australia 


CONTENTS 


Page  No. 


1.  INTRODUCTION  I 

1.1  Report  Outline  2 

2.  BASIC  RELIABILITY  FUNCTIONS  2 

2.1  Risk  Rate  KO  2 

2.2  Probability  of  Failure  P^t)  3 

2-3  Probability  of  Survival  P,(t)  3 

2.4  Probability  Density  for  the  Time  to  Failure  3 

2.5  Relationships  Between  the  Basic  Reliability  Functions  3 

2.6  The  Distribution  of  a  Structural  Characteristic  4 


3.  MATHEMATICAL  BASIS  FOR  A  RELIABILITY  MODEL  5 

3.1  A  General  Expression  for  the  Probability  of  Survival  S 

3-2  A  Basis  for  a  Fatigue  Model  6 

3.3  The  Effects  of  Inspections  8 

3.4  Generation  of  Derived  Functions  9 

3.4.1  Probability  of  Failure  9 

3.4.2  Probability  Desaity  for  tbe  Time  to  Failure  and  Risk  Rate  10 

3.4J  Failure  Density  for  Strength  11 

3.4.4  Probability  Density  for  Strength  12 

33  Some  General  Results  12 

3.5.1  The  Effect  of  a  Constant  Additional  Risk  12 

333  The  Effect  of  a  Ham  Dependent  Additional  Risk  13 

3  A3  Transformation  of  the  Random  Variables  14 

3.6  Summary  of  Steps  Required  to  Generate  a  Reliability  Model  for  Fatigue  16 


4.  THE  TWO-PARAMETER  MODEL  OF  PAYNE  ET  AL. 

4.1  Assumptions  and  Model  Equations 

4.1.1  Random  Variables 

4.1.2  Time  Zones  and  Physical  Processes 
4.13  Risk  Rate  Equations 

4.1.4  Subspace  Boundaries 
43  Transformations  of  tbe  Random  Variables 
43  An  Integral  Expression  for  P^t) 


n 

17 

17 

18 
18 
18 
20 
20 


HR 


wv 


S* - -  .  ..  .  >  „  V* 

Justification. 


By - 

Dlstrl but 1 oa/ 


Availability  Codas 


01st 


|Avall  and/or 
Special 


□  U 


4.4  Inspections 

23 

4.5  Generation  of  MO  •**  5-«k  Rate* 

24 

4.6  Probability  Density  for  Strength 

26 

4.7  Failure  Density  for  Strength 

27 

4.8  Approximate  Expressions  for  the  Risk  Rate 

28 

5.  THE  TWO-PARAMETER  MODEL  OF  HOOKE 

29 

5.1  Assumptions  and  Mode)  Equations 

29 

5.2  An  Integral  Expression  for  PJO 

29 

5.3  Calculation  of  />(/)  and  Risk  Rates 

30 

5.4  The  Effect  of  Fatigue  Life  limiting 

30 

5.5  Failures  Caused  by  Fatigue  Weakening 

31 

6.  THE  TWO-PARAMETER  MODEL  OF  FORD 

33 

6.1  Assumptions  and  Mode)  Equations 

33 

6.2  An  Integral  Expression  for  P*(t) 

35 

6.3  Calculation  of  p^r)  and  Risk  Rates 

35 

6.4  Transformed  Random  Variables 

36 

A5  Ford’s  One-Crack  Model 

39 

6.6  Inspections 

40 

6.7  Failure  Density  for  Crack  Length 

40 

6JI  Model  Extemkns  to  Iadade  Initial  Cracking 

40 

6A1  Initial  Age  as  a  Random  Variable 

41 

62.2  The  Initial  Crack  Model  of  Diamond  and  Payne 

41 

7.  A  THREE-PARAMETER  MODEL 

42 

7.1  Assumptions  and  Model  Equations 

42 

7.1.1  Bask  Random  Variables 

42 

7.1.2  Transformed  Random  Variables 

42 

7.1.3  Initial  Age 

43 

7.1.4  Model  Equations  and  Sabapnce  Boundaries 

43 

7.2  Integral  Expression  for  P«<0 

44 

73  Risk  Rate  Equations 

44 

7.4  Redaction  to  Simpler  Models 

44 

g.  CONCLUSIONS 

46 

REFERENCES 

NOTATION 

DISTRIBUTION 

DOCUMENT  CONTROL  DATA 

- •  w-i -Hir^-  riWi^laiiifri^fftii>irta#i»1‘  i'  *  -a^*-**,  ;*oi  «*iw5**‘  -.«  .«* 


1.  INTRODUCTION 

There  are  many  phenomena  (e.g.  fatigue,  corrosion),  which  progressively  degrade  the 
ability  of  a  structure  or  component  to  survive  the  effects  of  its  service  environment.  Ultimately 
this  degradation  results  in  the  complete  failure  of  the  structure  at  a  time  which  is  a  function  of 
the  environmental  and  degradation  histories.  Unfortunately,  the  uncertainties  in  both  the 
environmental  history  and  its  effect  on  the  degradation  process  make  a  prior  deterministic 
calculation  of  the  time  of  failure  impossible. 

Certain  classes  of  structures,  such  as  aircraft  and  automobiles,  are  manufactured  in  such  a 
way  that  for  a  given  environment,  the  degradation  processes  are  similar  between  members  of 
each  class.  It  is  then  possible  to  apply  statistical  methods  to  estimate  the  time  of  failure  based 
on  a  knowledge  of  the  mean  behaviour  of  the  class,  deduced  from  experiments  or  from  service 
failures.  Similar  statistical  methods  can  be  used  to  account  for  variations  in  the  environment 
and/or  the  effect  that  the  environment  has  on  the  degradation  process. 

A  reliability  model  is  one  such  statistical  method.  The  objective  of  the  model  is  the  compu¬ 
tation  of  reliability  functions  which  predict  the  behaviour  of  an  ensemble  of  structures  when 
subjected  to  sequences  of  applied  loads  (the  environment)  which  fall  within  given  statistical 
descriptions. 

There  are  two  approaches  which  may  be  used  to  generate  a  reliability  model  for  fatigue. 
The  first  approach  is  to  assume  that  strength  degradation  is  related  to  crack  length  which  is 
either  a  known  function  of  time  (for  a  given  load  application  rate)1-*  or  determined  by  a  dif¬ 
ferential  equation.9 ~li  Variations  in  structural  behaviour  are  accommodated  by  introducing 
random  variables  as  parameters  in  the  strength-crack  length-time  relationships. 

The  second  approach  is  to  assume  a  probabilistic  structure  from  the  start  and  describe  the 
evolution  of  the  probability  distributions  of  critical  descriptive  random  variables  using  appro¬ 
priate  statistical  techniques.  An  example  of  this  approach  is  the  Markov  Chain  model  described 
by  Bogdanoff  and  Kozin12  who  also  presented  a  critical  appraisal  of  the  two  approaches  for 
modelling  fatigue  crack  growth. 

The  models  developed  at  ARL  by  Payne  et  ai.,1  Hooke5-*  and  Ford9-11  have  used  the 
first  approach  which  has  the  advantage  that,  although  fatigue  is  by  no  means  a  deterministic 
process,  the  ensuing  model  has  the  appeal  of  having  a  more  apparent  physical  basis. 

In  parallel  with  the  progressive  development  of  the  reliability  models  at  ARL,  techniques 
for  their  numerical  evaluation  have  been  established  and  refined,  culminating  in  the  NERF 
(Numerical  Evaluation  of  Reliability  Functions)  computer  program  described  by  Mallinson 
and  Graham.13 

The  work  reported  herein  originated  during  the  development  of  the  NERF  computer  pro¬ 
gram  when  it  became  evident  that  existing  reports  describing  the  reliability  theory  presented 
neither  sufficiently  complete  analyses  nor  adequate  details  of  the  reliability  functions  to  provide 
a  basis  for  the  specification  of  the  numerical  procedures.  The  computer  program  was  designed 
for  the  Payne  et  ai.  models.  It  was  not  clear  if  the  relationship  between  those  models  and  the 
analyses  presented  by  Ford9- 11  was  sufficiently  strong  to  permit  application  of  the  numerical 
techniques  to  Ford’s  models.  Moreover,  despite  similar  assumptions,  the  models  of  Payne  and 
Hooke  were  not  in  complete  agreement*  regarding  the  expressions  for  some  of  the  reliability 
functions. 

As  this  theoretical  study  progressed,  it  became  evident  that  a  unified  method  for  the  deriva¬ 
tion  of  reliability  models  could  successfully  yield  all  the  reliability  functions  developed  at  ARL 
and  similar  functions  used  elsewhere.14-15  Differences  between  the  various  models  could  be 
interpreted  and  the  method  fulfilled  the  original  objective  of  generating  a  detailed  specification 
for  the  NERF  computer  program.13 

It  should  be  stressed  here  that  the  investigation  and  this  report  considers  only  the  mathe¬ 
matics  behind  the  derivation  of  reliability  models.  It  does  not  consider  the  problem  of  obtaining 
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suitable  input  data  for  the  model.  Hie  various  functions  and  probability  density  functions  used 
by  the  reliability  analysis  are  assumed  to  be  known.  The  estimation  of  the  various  parameters 
required  by  any  model  for  fatigue  is  non-trivial  and  requires  very  careful  analysis  of  the  available 
data,16-17  and  a  complete  investigation  of  the  accuracy  of  such  estimations  and  the  effect  that 
any  errors  may  have  on  the  ability  of  the  reliability  model  to  predict  a  “real  life”  situation  was 
beyond  the  scope  of  this  investigation. 

1.1  Report  Outline 

Following  the  definition  of  the  basic  functions  evaluated  by  a  reliability  model  (Chapter  2), 
the  mathematical  basis  of  a  unified  method  for  the  generation  of  suitable  mathematical  expres¬ 
sions  for  these  functions  is  described  (Chapter  3).  The  analysis  indicates  how,  for  any  reliability 
model,  the  effects  of  inspections  can  be  included  and  density  functions  for  any  desired  structural 
characteristic  obtained.  Without  making  any  assumptions  regarding  the  precise  form  of  the  model, 
general  results  such  as  the  effect  of  additional  but  independent  hazards  to  structural  survival  are 
deduced.  In  many  cases,  it  may  be  advantageous  to  make  a  transformation  from  the  original 
random  variables  used  in  the  model  to  new  variables  to  obtain  forms  of  the  reliability  functions 
which  are  easier  either  to  evaluate  or  to  compare  with  those  associated  with  a  different  model. 
The  mathematics  associated  with  effecting  such  transformations  are  described.  Chapter  3  ends 
with  a  summary  of  the  steps  required  to  generate  a  reliability  model  from  the  initial  assumptions. 

The  remainder  of  the  report  demonstrates  the  application  of  the  method  described  in  Chapter 
3  to  derive  specific  reliability  models.  Because  of  the  relevance  to  the  documentation  of  the  NERF 
computer  program,  the  two-parameter  model  developed  by  Payne  et  al.l~A  is  considered  first 
(Chapter  4)  in  some  detail.  The  original  random  variables  used  by  Payne  et  al.  are  transformed 
to  new  variables  representing  “age”  and  “virgin  strength”.  The  resulting  model  which  is  a  little 
more  detailed  than  the  Payne  model  verifies  the  functions  reported  by  Payne  and  Graham.4 
The  relevance  of  approximate  expressions  derived  by  Diamond  and  Payne3  is  considered. 

The  closely  related  model  developed  by  Hooke5-8  is  studied  in  Chapter  5.  In  terms  of  the 
transformed  random  variables  the  differences  between  the  models  of  Payne  and  Hooke  are 
elucidated.  Hooke’s  proposal  to  isolate  failures  by  fatigue  weakening  is  examined. 

Ford’s  models’ were  developed  using  methods  different  from  those  used  by  Payne  and 
Hooke.  Expressions  identical  to  those  presented  by  Ford9,11  are  generated  by  the  method 
described  here  (Chapter  6).  Ford’s  model  is  extended  to  include  structures  that  are  initially 
cracked  and  the  similarity  between  this  extended  model  and  the  initial  crack  model  proposed 
by  Diamond  and  Payne3  is  investigated. 

A  three-parameter  model  is  developed  in  Chapter  7.  This  model  encompasses  all  the  models 
described  in  the  previous  sections  of  the  report  as  more  simple  cases.  A  tabulation  of  each  model 
in  terms  of  the  three-parameter  model  is  presented. 


2.  BASIC  RELIABILITY  FUNCTIONS 

The  objective  of  a  reliability  model  is  to  evaluate  “reliability  functions"  which  represent  the 
aggregate  behaviour  of  a  population  of  structures  which  fall  within  the  scope  of  a  given  statistical 
representation.  These  basic  reliability  functions  are  defined  below. 

2.1  Risk  Rate  KO 

The  overall  effect  of  the  environment  and  degradation  processes  is  the  depletion  of  the 
population  at  a  time  dependent  rate.  Defining  the  removal  of  a  structure  from  the  population 
as  a  “failure”,  this  rate  can  be  represented  as  a  risk  rate  defined  in  words  by, 

r\t)  =  risk  rate 

__  No.  of  failures/(unit  time) _  ^ 

~  No.  of  structures  in  population  at  time  t 

The  superscript  denotes  a  definition  based  on  a  population  with  a  finite  number  of  members.  As 
this  number  increases,  the  risk  rate  approaches  the  continuous  definition, 

r(t)  =  risk  rate  —  Fraction  of  remaining  population  failing/(unit  time),  at  time  /.  (2.2) 
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22  Probability  of  Failure  P,(t) 

Let  the  random  variable  F  denote  the  time  of  failure  of  a  structure  in  the  population.  At 
time  t,  the  probability  of  failure  is  the  probability  that  F  is  less  than  t.  This  probability  can  be 
written  as  a  cumulative  distribution  Pf(t)  where. 


P%t)  =  P(F  <  r)  = 


Number  of  failures  up  to  time  t 
Original  number  of  structures 


or  for  a  large  population. 


(2.3) 


Pi(t)  =  Fraction  of  original  population  that  has  failed  before  time  t.  (2.4) 

Pf(t)  ranges  from  0  when  t  =  0  to  1  as  /  becomes  infinite. 


12  Probability  of  Survival  J*g(t) 

At  time  t,  the  probability  of  survival  is  the  probability  that  F  is  greater  than  f.  For  a  large 
population, 

P^t)  =  P{ F  >  t)  =  Fraction  of  original  population  remaining  at  time  f.  (2.5) 

The  total  population  must,  at  any  time,  be  composed  of  two  mutually  exclusive  sets  of  structures, 
those  that  have  failed  and  those  that  have  survived.  Obviously, 

W+JVfl-1.  (2.6) 


2.4  Probability  Density  for  the  Time  to  Failure  #*(t) 

The  fraction  of  the  population  with  times  of  failure  in  the  interval  ( r.t+dt )  is  given  by 
Pw{t)dt,  where  pjj)  is  the  probability  density  function  for  the  time  of  failure. 


2-5  Relationships  Between  the  Basic  Reliability  Functions 

From  the  definitions  given  above,  relationships  between  the  reliability  functions  can  be 
derived.  The  probability  distribution  of  the  time  to  failure  is  related  to  p^t)  by 


P* 0  = 


from  which  it  follows  that 


P»(  0  = 


dPM  -dP#) 
dt  dt  ' 


For  a  finite  time  interval,  At,  the  number  of  failures  is 

p/.t) .  At  x  (the  number  of  structures  in  the  population). 

It  follows  that, 

_  MO 


no  * 


and  in  the  limit  as  n  -*■  co, 


dPdLOl. 

dt  / 


Th 


3 


or 


Equation  (2.10)  can  be  integrated  to  yield, 


w= 


(2.10) 


(2.11) 


2.6  The  Distribution  of  a  Structural  Characteristic 

Let  the  random  variable  Z  denote  a  particular  structural  characteristic.  The  fraction  of  the 
population  with  z<Z<z+dz  failing  in  the  interval  (tj+dt)  is  given  by  p,T(t^:)dt .  dz  where 
Pr.z(‘>z)  *s  the  joint  density  function  for  F  and  Z.  If  Pr|Z(r,z)  is  integrated  over  all  values  of  z,  the 
total  fraction  of  the  population  failing  in  the  interval  (t,t+dt)  is  obtained,  i.e., 

f°° 

Pf(0  =  Pr.z«,z>fc.  (2.12) 


The  probability  density  of  the  time  to  failure  is  the  marginal  density  for  F. 

At  a  particular  value  of  time,  it  can  be  of  considerable  interest  to  know  the  distribution 
of  Z  among  the  failing  structures.  This  information  is  obtained  by  evaluating  the  conditional 
density  for  Z  given  failure  at  time  f,  pz(z|F  =  f),  written  here  as  pz(z|/)  where. 


Pt,z{t,z)dz 


— oo 


Pr.Z(t,z) 

P/t) 


(2-13) 


/’oo 

Conversely,  if  p^t)  can  be  expressed  in  the  form  I  f{z)dz,  then, 

J  -oo 


PZ(z|0  =/(z)/Pr(0- 


(2.14) 


The  conditional  density  for  Z  given  failure  at  time  t  will  be  referred  to  here  as  the  failure  density 
of  the  characteristic. 

At  time  t,  the  fraction  of  the  population  with  z  <  Z  <  z+dz  that  has  survived  is  given  by 
*oo 

PftZ(t,z)dt  |  dz.  The  total  fraction  surviving  to  time  t  is  given  by  the  marginal  probability 


If 


of  survival. 


Pi.zkt,z)dtdz  =  PJt). 


The  conditional  density  for  Z  given  survival  to  time  t  is  pj(z|F  >  /)  where, 


Pa(z|F  >  0 


W) 


(2.15) 


(2.16) 
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Conversely,  if  /*(/)  can  be  written  in  the  form  I  g(z)dz, 

J  -oo 

pM*  >  0  =  &)/?£)■  (2.17) 

This  density  function  will  be  referred  to  here  as  the  density  function  for  the  characteristic. 


3.  MATHEMATICAL  BASIS  FOR  A  RELIABILITY  MODEL 

The  derivation  of  a  reliability  model  for  fatigue  can  be  performed  by  several  methods. 
Payne  et  al.l~ 4  derive  the  aggregate  risk  rates  from  first  principles.  Ford9-11  derives  an  equation 
of  continuity  for  the  probability  density  of  crack  length  and  integrates  this  equation  to  generate 
the  reliability  functions.  Hooke*  “•  derives  both  the  probability  of  survival  and  the  average 
risk  rates  from  first  principles. 

The  approach  used  here  is  to  derive  the  probability  of  survival  directly  and  then  use  the 
relationships  established  in  the  previous  section  to  generate  the  remaining  reliability  functions. 
This  approach  resembles  that  of  Hooke,*-*  but  is  more  general. 


3.1  A  General  Expression  lor  the  Probability  of  Survival 

The  derivation  of  an  expression  for  the  probability  of  survival  starts  with  the  following 
assumptions. 

(i)  The  ]  filiation  of  structures  can  be  characterised  by  a  set  of  M  random  variables, 

X2, . . .,  XM  with  joint  probability  density  function  />*(x)  where,  to  reduce  subse¬ 
quent  notation,  X  denotes  the  vector  of  random  variables  Xi,  X2, . . .,  XM  and  x 
denotes  the  vector  of  sample  values  xu  x2, . . .,  xM. 

(ii)  A  given  structure  has  a  particular  combination  of  values  of  the  random  variables. 
These  values  remain  constant  for  the  life  of  the  structure. 

(iii)  Each  structure  passes  through  a  sequence  of  time  zones  during  each  of  which  a  distinct 
physical  process  can  be  represented  by  an  appropriate  risk  rate  that  is  a  function  of  the 
random  variables  and  time. 

(iv)  The  number  and  ordering  of  the  time  zones  is  the  same  for  every  member  of  the 
population  and  each  structure  can  fail  as  a  result  of  only  one  of  the  physical  processes. 
(The  models  described  herein  are  strictly  one-crack  models.) 

Expressed  formally,  for  given  x  there  exists  a  sequence  of  times,  f„(x),  /,(x), . . .,  r*(x)  such 
that  for  /*_!</<  tk;  (k  <  K), 

r(t)  =  r»(x,f)  *  *■*(*„  x2, . . .,  xM,  t).  (3.1) 

Conversely,  for  a  given  value  of  time,  there  exists  a  set  of  subspaces  {Dk}  of  x  space.  In  the 
subspaces  Dk 

fit)  =  rk(x,t). 
x 

The  subspaces  are  disjoint,  i.e.  Dif)D,  =  0  and  U  Dk  =  whole  of  x  space.  As  time  progresses 
a  structure  with  given  x  moves  sequentially  through  the  subspaces. 

Consider  those  structures  belonging  to  the  M  dimensional  hypervotume 
(*1,*i+<(*i)x(*i»*a-l'«fcri)x.  •  .x(**,x* +<&*). 


The  fraction  of  the  original  population  that  has  these  values  of  the  random  variables  is 

m 

Px(x)dx  (where  dx  denotes  FI  <**<)•  As  time  proceeds,  this  fraction  of  the  population  is  depleted 

lal 

by  the  successive  risk  rates,  rt (x^);  1  <  k  <  K.  Because,  in  this  volume  of  x  space  the  random 
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variables  are  essentially  constant,  equation  (2.13)  can  be  used  to  yield  the  probability  of  survival 
and  hence  the  fraction  of  this  section  of  the  population  remaining  at  time  t.  For  <  t  <  tk, 
and  denoting  the  local  probability  of  survival  by  the  result  of  this  integration  is. 


-  Z  rjx,t)di-  I  r^t)dt 

'“J  tj-i  J  4-i 

The  surviving  fraction  of  this  section  of  the  original  population  is, 

P^x,t)p^x)dx  s  dFk(x,t),  say. 


P%x,t)  =  exp 


(3.2) 


(33) 


The  total  fraction  of  the  original  population  that  has  survived  to  time  t  can  be  obtained  by 
integrating  (3.3)  over  the  whole  of  x  space,  remembering  that  there  are  K  distinct  subspaces  to 
integrate  over. 


r/x,t)dt—  f  rk(t)dt 

J  4-i 


(3-4) 


(Note  that  the  integration  domain  for  each  Dk  is  time-dependent.)  This  expression  for  the 
probability  of  survival  is  general  and  can  be  used  to  derive  a  reliability  model  for  any  process 
which  satisfies  the  assumptions  made  above  and  for  which  suitable  local  (i.e.  in  terms  of  the 
random  variables)  risk  rates  can  be  found. 


3  J.  A  Basis  for  a  Fatigue  Model 

Hooke7  and  Saunders17  identify  three  distinct  time  zones  before  failure.  In  this  discussion, 
two  are  identified  and  are  considered  to  be  sufficient  to  provide  an  adequate  model  within  the 
scope  of  available  fatigue  data.  The  resulting  total  of  three  time  zones  (including  failed  structures) 
is  defined  below. 

(i)  Dk  :  Uncracked  structures 

From  initial  time  (/„)  until  the  instant  of  crack  initiation,  the  structure  is  assumed  to 
be  uncracked  and  its  strength  the  same  as  when  it  was  new,  or  in  its  virgin  state.  There 
may  be  some  risk  of  failure  due  to  service  loads  being  greater  than  the  virgin  strength. 

(ii)  Z)2:  Cracked  structures 

At  time  tu  a  fatigue  crack  is  assumed  to  initiate,  following  which  the  structure  is 
weakened  and  suffers  an  increased  risk  that  applied  loads  will  exceed  its  strength.  For 
the  purpose  of  this  analysis,  crack  initiation  is  assumed  to  correspond  to  the  time  at 
which  strength  starts  to  reduce.  The  cracked  but  unweakened  phase  identified  by  Hooke 
is  considered  to  be  indistinguishable,  from  a  risk  point  of  view,  from  the  uncracked 
phase. 

During  the  second  time  zone,  the  structure’s  strength  continues  to  reduce  as  the 
fatigue  crack  grows. 

(iii)  Ds:  Failed  structures 

A  frequently  adopted  assumption  of  a  reliability  model  is  that  there  is  a  finite  fatigue 
life  for  any  structure.  Thus,  in  addition  to  removal  from  the  population  as  a  result  of 
an  applied  load  bring  greater  than  the  strength  of  the  structure,  the  population  is 
depleted  by  structures  reaching  their  fatigue  life.  This  event  is  often  associated  with 
extremely  rapid  crack  growth  at  t  «■  t2  (Ford)  or  with  a  sufficient  reduction  in  the 
strength  of  the  structure  to  make  it  fail  under  an  applied  constant  load  (Payne). 

Note  that  despite  the  name  given  to  the  time  zone,  D3  does  not  contain  all  the 
structures  that  have  Med;  only  those  that  have  reached  the  fatigue  time  limit. 
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Clearly,  the  total  fraction  of  the  original  population  remaining  at  any  time  involves  integration 
over  Dx  and  D2  only. 

To  generate  a  complete  model,  an  expression  for  the  local  risk  rate  in  each  time  zone  must 
be  derived.  In  general,  it  is  in  this  area  that  the  differences  between  various  fatigue  models 
emerge.  There  is,  however,  a  certain  amount  of  the  analysis  that  can  be  made  before  reaching 
the  model  dependent  stage. 

In  the  case  of  fatigue,  the  risk  rate  is  assumed  to  be  the  consequence  of  an  applied  fluctuating 
load  sequence  which  has  two  major  effects. 

(i)  Damage  assumption 

Each  load  application  produces  an  increment  of  damage  which,  if  the  structure  is 
uncracked,  produces  an  inevitable  approach  to  crack  initiation.  Following  initiation, 
the  damage  becomes  an  increment  in  crack  length  and  a  subsequent  reduction  in 
strength. 

(ii)  Contribution  to  risk 

Each  load  application  has  the  potential  for  producing  structural  failure  if  the  load 
exceeds  the  current  strength  of  the  structure. 


These  two  effects  can  be  combined  to  generate  a  risk  function  in  the  following  manner. 
The  damage  assumption,  and  the  fact  that,  in  terms  of  the  fatigue  life  of  a  structure,  the  applied 
loads  occur  with  sufficient  frequency  that  variations  in  their  timing  may  be  ignored,  permits  a 
direct  link  between  time  and  load  applications  to  be  made  via  an  average  load  application 
rate,  /,. 

Let  R  denote  the  current  strength  of  a  structure  and  let  L  denote  load  magnitude.  From  a 
given  load  sequence  it  is  possible  to  regard  I  as  a  random  variable  and  derive  a  probability 
density  pJL)  for  the  applied  load.  The  probability  distribution  of  applied  loads  is. 


PdL)  = 


Pt(L')dL'. 


(3-5) 


Assuming  that  the  applied  load  sequence  is  statistically  invariant  with  respect  to  time,  the 
probability  that  an  applied  load  exceeds  the  current  strength  of  a  structure  is 

P(L  >  R)  =  1  ~Pl(R)  =  Ph(R).  (3.6) 


The  instantaneous  risk  rate  is  given  by 

'2(0  =  *«.(*)•/,•  (3-7) 

The  damage  assumption  is  included  by  postulating  that  the  strength  of  a  cracked  structure 
is  a  function  of  crack  length  which  is,  in  turn,  a  continuous  function  of  the  number  of  load 
applications.  This  can  be  represented  formally  by, 

R  =  tf(a(0)  (3.8) 

where  a  denotes  crack  length. 

For  an  uncracked  structure  the  strength  is  assumed  to  be  constant,*  Ro  say,  so  that, 

r,0)  =  Pt(Ro)  ■  /r-  (3-9) 

These  expressions  for  risk  pertain,  of  course,  to  either  a  particular  structure  or  to  a  uniform 
population  of  structures.  Variations  between  structures  can  be  accounted  for  by  introducing 
into  equations  (3.7H3.9)  parameters  which  become  the  random  variables  of  equation  (3.1).  It 
is  possible  to  introduce  such  parameters  in  a  completely  general  way  thereby  accounting  for 
variations  in  the  applied  load  sequence,  crack  growth  and  strength  decay.  The  models  that  have 
been  used  to  date,  and  those  described  herein,  consider  statistical  variations  in  the  last  two 
mechanisms  only.  Random  variables  need  only  be  introduced  into  the  expressions  for  R,  viz., 

Ro  -  Ro(a),  (310) 


*  There  is  evidence,  however,  that  strength  may  reduce  before  crack  initiation.1* 
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R  =  R(xj). 


(3.U) 


(Note  that  in  general  R0  #  R(x,t0).) 

To  simplify  subsequent  notation,  the  following  equivalent  expressions  will  be  used  with  the 
simplest  being  preferred  where  possible  without  ambiguity. 

'i  =  ri(x)  =  PM*)) .  /„  (3.12) 

r2  =  r2(x,t)  =  PL(R(x,t)) .  /r.  (3.13) 

The  boundaries  between  the  three  subspaces  are  defined  by  the  equations, 

/  =  f,(x)  for  DJD2,  (3.14) 

t  =  t2(x)  for  D2ID2,  (3.15) 

where  DJD2  denotes  the  boundary  between  Z),  and  D2,  etc. 

The  analysis  leading  to  equation  (3.5)  may  now  be  followed  through.  Equation  (3.2)  leads  to, 

Pa(x,t)  =  exp{— r,(x) .  r},  (3.16) 


P i(x,0  =  exp  — r 


r*  \ 

ri(x)  •  h-  r2(x,t)dt  J, 
J  h  ) 


so  that  the  expression  for  the  probability  of  survival  becomes 


Pa(‘)=  j  exp  {-r1t}px(x)dx+  expl-r^j-  r2(x,l)di  Lx(x>/x.  (3.18) 

J  D  i  D2  [  Jh  j 

A  final  simplification  to  the  notation  can  be  made  by  defining  a  loss  factor 

•  •  •*  Oit 
(Pl(*,0  for  t  < 


H{x,t)~ 


P\(x,t)  for  ?!  <  /  <  t2. 


which  upon  replacement  in  (3.18)  yields, 


P*(t)  = 


H(x,t)px(x)dx+  H(x, 

J  Dt  J  D2 


t)px(x)dx. 


Equations  (3.18)  or  (3.19)  and  (3.20)  define  a  probability  of  survival  from  which  all  other 
reliability  functions  can  be  generated  for  a  fatigue  model. 


3 13  The  Effects  of  Inspections 

An  important  requirement  of  a  reliability  model  is  the  ability  to  account  for  inspections  and 
subsequently  provide  a  basis  for  the  selection  of  inspection  strategies. 

An  inspection  has  two  effects  on  the  population  of  structures: 

(i)  structures  that  are  detected  to  be  deficient  according  to  appropriate  criteria  are  removed ; 

(ii)  the  structures  removed  during  an  inspection  may  be  replaced. 

The  lint  effect  is  relatively  easy  to  incorporate  into  the  reliability  model.  Following  an  inspection 
at  tii  ^y.  Pxb 0  can  be  replaced  by  p2(x)  where, 


t  Note  the  notational  difference  between  the  loss  factor  defined  here,  and  H,  used  elsewhere 


to  denote  the  integrated  Hazard 
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P&x)  =  Px(*)  •  soffit)  (3.21) 

and  S(x,ft,)  is  a  function  representing  the  removal  of  structures  during  an  inspection.  In  general, 
this  function  can  represent  both  the  application  of  the  inspection  testing  procedures  and  the 
probability  that  those  procedures  may  fail  to  detect  a  deficient  structure. 

The  probability  of  survival  following  the  inspection  is, 

P,(t)  =  H(x,t)pHx)dx+  H(x,t)px(x)dx  (3.22) 

J  D,  J  D2 

and  the  total  fraction  of  the  population  removed  during  the  inspection,  called  here  the  probability 
of  detection  at  ft„  is 

PWh)  =  Pi(th)~Pi(tii)  (3.23) 

where  /»)(«,)  and  Pj(fi,)  denote  the  probability  of  survival  evaluated  by  equations  (3.20)  and 
(3.22)  respectively. 

The  second  effect  is  more  difficult  to  incorporate.  In  general,  the  replacement  structures 
come  from  a  population  which  has  different  statistical  properties  from  the  original  population. 
At  best  the  replacements  may  be  selected  at  random  from  a  population  similar  to  the  original 
one.  However  the  replacements  may  be  subjected  to  repair  or  inspection  procedures  which 
would  alter  the  distributions  of  the  random  variables.  Accordingly,  the  probability  of  survival 
following  the  first  inspection  can  be  written  as, 

p,0)  =  P°(t)+P\t-til)  .  PUti,)  (3.24) 

where  P°(t)  represents  the  probability  of  survival  for  the  original  population  (equation  (3.22)) 
and  Px{t— tit)  represents  the  probability  of  survival  for  the  population  of  replacement  structures 
introduced  at  time  ft',. 

Following  a  second  inspection  at  ft)  say,  the  reliability  model  must  account  for  three 
populations;  the  original  population  where  pj(x)  is  now, 

p5(x)  =  PxOO‘S(x,«1)S(x,rt1),  (3.25) 

the  ft,  replacement  population,  similarly  reduced  by  S(x,ft2),  and  a  new  replacement  population 
of  size  Pd„(ft2).  As  further  inspections  occur,  the  number  of  populations  included  in  the  model 
continues  to  increase. 

Fortunately,  useful  calculations  can  be  made  without  involving  the  full  complexity  described 
above.  For  example,  in  the  models  described  by  Payne  et  al.,  the  inspections  are  assumed  to  be 
perfect  so  that  the  inspection  operator  S(x,t)  can  be  replaced  by  modified  integration  limits  (or 
subspace  boundaries).  Replacement  structures  are  assumed  to  be  repaired  in  such  a  way  that  they 
are  no  longer  susceptible  to  a  risk  of  failure  (Payne2)  or  the  calculations  may  proceed  without 
replacements  (Hooke*).  In  either  case  only  one  population  need  be  considered  in  the  analysis. 

In  this  report,  the  inspection  analysis  developed  by  Payne  et  al.  only  will  be  discussed.  A 
more  detailed  treatment  of  the  general  inspection  analysis  has  been  presented  by  Ford.10 


3 A  Generation  of  Derived  Fractions 

Given  the  expression  for  the  probability  of  survival,  those  for  the  remaining  basic  reliability 
functions  can  be  generated  by  straightforward  applications  of  the  relationships  (2.6)  to  (2.11). 
The  general  steps  involved  are  outlined  in  this  section  for  each  of  the  basic  functions  together 
with  the  probability  densities  for  strength  and  failing  load  which  are  defined  below  and  have 
particular  relevance  to  a  fatigue  model. 


W-i-JVM. 


3.4.1  Probability  of  Failure 
Using  (2.6), 
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(3.26) 


and  the  probability  of  failure  can,  in  principle  be  obtained  in  this  way.  In  practice,  however, 
numerical  problems  arise  since  a  high  degree  of  accuracy  in  the  determination  of  P^t)  is 
required  (particularly  for  low  values  of  t  when  Pt(t)  «  1),  to  ensure  that  a  realistic  estimate  of 
Pf(t)  is  made.  It  is  preferable  to  calculate  Pf(t)  directly  using  an  expression  analogous  to  (3.20). 

Remembering  that  Pf(x,t)  =  1  —  /»(x,f)  so  that,  in  particular,  P,  =  1,  the  expression  for 

Wis, 

PM=  I  (1  -H(x,t))px{x)dx+  I  (1  —H(x,t))px(x)dx+  I  px(x)dx.  (3.27) 

J  D  i  J  D  2  J  D3 

Despite  its  greater  complexity,  this  expression  is  in  practice  numerically  easier  to  evaluate  than 
(3.20).  The  computer  program  NERF  evaluates  P^t)  this  way  and  derives  Pt(t)  using  (3.26). 


3.4.2  Probability  Density  for  the  Time  to  Failure  and  Risk  Rate 

The  expression  for  the  density  function  for  the  time  to  failure  arises  from  the  second  relation¬ 
ship  in  equation  (2.8),  viz., 


M0  = 


~dPn(  0 
dt  ' 


(3.28) 


Before  the  details  of  the  differentiation  are  discussed,  it  is  worth  considering  the  physical 
meaning  of  the  integral  terms  in  equation  (3.20).  Each  term  is  an  expression  for  the  total  fraction 
of  the  population  contained  within  a  particular  subspace  and  the  time  differential  yields  the  rate 
of  change  of  that  fraction.  With  the  progression  of  time,  the  number  of  structures  in  a  subspace 
changes  for  two  reasons.  First,  the  risk  incorporated  in  the  loss  factor  implies  a  global  removal 
over  the  subspace.  Second,  the  boundaries  of  the  subspaces  are  time-dependent,  meaning  that 
structures  leave  or  enter  the  subspace  as  the  boundaries  sweep  the  x  space. 

For  each  subspace,  equation  (3.28)  can  be  expected  to  yield  three  terms;  the  first  representing 
losses  resulting  from  the  applied  risk  rate,  the  second  and  third  representing  the  flow  of  structures 
through  the  boundaries  of  the  subspace.  Generally,  the  second  and  third  terms  in  a  subspace  will 
exactly  cancel  similar  terms  from  other  subspaces.  For  example,  structures  crossing  the  initiation 
boundary  leave  Dv  and  enter  D2  with  no  nett  loss  or  gain  during  the  transition.  The  exception  is 
the  term  in  D2  which  represents  removal  into  the  subspace  D3  of  failed  structures.  This  term 
constitutes  a  direct  contribution  to  p/f)  and  subsequently  to  r(f).  It  is  the  source  of  the  “risk  of 
fatigue  fracture’*  introduced  by  Payne  and  Grandage1  and  later  questioned  by  Hooke.7 

To  simplify  subsequent  analysis,  but  nevertheless  demonstrate  the  essential  mathematics, 
consider  a  model  which  has  two  independent  random  variables  with  density  functions  PXl(x i) 
and  PXl(x2).  The  subspace  boundaries  are  assumed  to  be  such  that  the  term  for  D2  in  (3.20)  can 
be  written  in  the  form, 

f*2,2  f  *1,2(0 

r2(0=  H(xlrxi,t)pXl(xl)pXl(x2)dxldx2  (3.29) 

J  *2,1  J  *i,i(0 

where  Xj.j  and  x21  are  solutions  for  equation  (3.14),  xl>2  and  x2>2  are  solutions  for  equation 
(3.1S)  and  only  Xj.,  and  x(>2  are  functions  of  time.  The  assumption  that  the  limits  of  the  inner¬ 
most  integration  are  the  only  limits  that  are  time-dependent  is  certainly  valid  for  the  models 
considered  in  this  report.  The  time  differential  of  this  term  can  be  written  as, 

<fTj(0  d 

^</l(*2.0)ftl1(*2)^2  (3.30) 

J  *2,1 

where 

f*  i,s(0 

A(*  2.0=  I  H(x  i^2,<)Px,(*i)^*i.  (331) 

J  *i.i(0 
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<u1_ 

*“J*  ,40)5 


tf(*i.*2.0Pn,(*i)*ri +dx-^~H{xx  ,i^t)PxSxx.i)- 


*i.i(0 


*(*i,iA>fW*i,i) 


-H(xltX2,t)  =  -r,(x1pc1,f)/r(x„x2,r)  (for  4-i  <  l  <  fj  (3.33) 

so  that 

rfT,  f*2-2  C*** 

"  J  *2.1  J  *„(!) 

f Xj-2  rfx,  ,(f) 

J  *2.1 

f Xj-2  dx,  ,(/) 

—  I  — ^  (3.34) 

J  *2.1 

By  generating  a  similar  expression  for  dTJdl,  it  can  be  readily  verified  that  the  third  term  in 
(3.34)  cancels  with  an  identical  term  representing  the  flux  of  structures  from  Dx  across  the  crack 
initiation  boundary. 

The  expression  for  Pf(t)  becomes, 

f*z.i  C *i.i 

=  ri(xltxt)H(xl^2,t)pXl(xi)p%1(x1)dxxdx2+ 

J  *2.0  J  *1,0 
f* 2.2  f*  1.2 

+  ri(xl^1,t)H(xl,x1,t)p1.Xxl)pXl(x2)dxldx2- 

J  x2,l  J  *1,1 

-j*  “p^(*i,ar*2.0Px,(*i,2)Px,<*2>*fj-  (3-35) 

Noting  that  r(t)  =  p/t)P^t),  then  r(r)  can  be  written  as  the  sum  of  three  risks, 

r(t)  =  rJj)+rt(t)+rAt\  (3.36) 

where  ry(f),  r,(t)  and  r/t)  correspond  to  the  three  terms  in  equation  (3.35)  and  are  named  here 
the  virgin  risk,  the  risk  of  static  failure  by  fatigue  and  the  risk  of  fatigue  life  exhaustion.  (The 
second  name  corresponds  with  that  used  by  Payne  el  al. ;  the  third  corresponds  with  their  “risk 
of  fatigue  fracture”.) 


ri(xlSi,t)H(xl,x1,t)p1Li(x1)pXl(x2)dxldx2- 


3.43  Failure  Density  of  Strength 

At  any  time  it  can  be  of  considerable  interest  to  know  the  magnitude  of  the  loads  most  likely 
to  cause  failure.  This  information  can  be  obtained  by  treating  the  strength,  A,  as  a  random 
variable  and,  following  the  analysis  of  Section2.6,  compute  the  failure  density  of  strength  (also 
called  probability  density  of  the  failing  load  by  Diamond  and  Payne*).  If  the  integral  expression 
for  Pf(t)  can  be  transformed  to  the  form 

MO- J*JaA)«  (3.37) 
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using  (3.1 1),  then  the  required  density  is. 


3.4.4  Probability  Density  for  Strength 

Similarly,  the  conditional  density  for  R  given  survival  to  time  t  can  be  evaluated  by  trans¬ 
forming  the  integral  expression  for  Pt(t)  to  the  form, 


P„0)  =  I  g(R)dR, 


^|F>0  =  ~-  (3-40) 

Examples  of  the  transformations  required  to  generate  the  densities  defined  by  equations  (3.38) 
and  (3.40)  are  given  in  Section  4  for  the  Payne  model. 


3.5  Some  General  Results 

Before  summarising  the  steps  required  to  generate  a  reliability  model  for  fatigue,  it  is  worth¬ 
while  considering  some  of  the  general  results  that  can  be  deduced  from  the  model  equations  as 
they  stand,  and  therefore  apply  to  any  of  the  models  derived  herein. 


3.5.1  The  Effect  of  a  Constant  Additional  Risk 
Let 

r{(0  =  rl(x)+rc,  (3.41) 

'i(  0  =  r2(x,t)+rc,  (3.42) 

where  r„  is  a  constant  risk  rate  probably  resulting  from  a  process  different  from  fatigue,  e.g.  the 
hijack  risk  identified  by  Ford.9  Equation  (3.2)  yields, 

P;‘(x,0  =  exp{—  [r,(x)+rjf}  =  Pj(x,f)exp{-rc»}  (3.43) 


Pg\x,t)  =  exp  -r,(x; 


r2(x,f)d<-rcf,-rc(<-f,) 


=  ?5(x,t)exp{-ret}. 

From  (3.43)  and  (3.44)  it  is  clear  that, 

/>#)  =  *xp{— r,t}Pg(»). 

Now, 


so  that 


jW-  -'«J>s(0+exp{-r,f}^/»i(0, 


PiiO  =  r.exp{  -r,/}P<(i)+exp{ -ref  }Mt) 
=  exp{ 
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r'{t)  =  r.+Mt)IP* 0 

=  rc+r(t). 

(3.47) 

If 

'w' 

II 

®  8 

1 

then 

Pft)  =  Jjs'(W 

where 

*'(*)  =  exp{-rc/}g(J?). 

(3.48) 

If 

1 

8 _ © 

II 

* 

then 

/*/)  =  Jj/'WdK 

where 

f\R)  =  exp{-re/}U(*)r,+/(Jl)]. 

(3.49) 

It  follows  that. 

P^R\¥>t)  =  ^=p^R\¥>,) 

(3.50) 

and 

rcg{R)+f(R) 

rcP4t)+M»' 

(3.51) 

3A2  The  Effect  of  a  Tine  Dependent  Additional  Utah 

Let 

ri(0  =  r,(x)+r.(t) 

(3.52) 

and 

ri(0  =  r2(x,J)+r,(t) 

(3.53) 

where  r,(<)  is  an  additional  risk  which  although  time-dependent,  is  not  a  function  of  the  random 
variables.  Equation  (3.2)  yields, 

Pi\x,t)  =  exp  —r  ,(*)/—  j  r,(t)dt 

Jo 

PiUx,  f)exp 


t 

rtm 

0 


(3.54) 


and 


Pi2(x,t)  =  exp 


from  which. 


Now, 


so  that 


and 


f'  f*i  f' 

-'■iW'i-  r2(x,t)dt-  I  r.(0</f-  r.(t)dt 

J  'i  JO  J  *1 


=  Pj(x,t )cxp[  -  I  r.(l)rfrj 


J*s(0  =  exp 


=  -r.(/)/>'(/)+exp{  -  |  rt(t)dt 


-j\(/>*J/>i('). 

ft 


Pift)  =  exp 


fUl'- 


MP'M+MO 


(3.55) 


(3.56) 


(3.57) 


r'(t)  =  r.(t)+r(r).  (3.58) 

The  results  for  p^(/t|F  >  l)  and  pi(/J|/)  are  identical  with  equations  (3.50)  and  (3.51)  (with  rc 
replaced  by  rj  for  the  case  of  constant  additional  risk. 


3dJ  Transformation  of  the  Random  Variables 


In  some  models,  it  may  be  convenient,  particularly  from  the  viewpoint  of  their  numerical 
evaluation,  to  recast  the  reliability  functions  in  terms  of  a  new  set  of  random  variables 

{Y„  Yj, . . .,  YM}  which  are  related  to  the  original  set  {Xt,  X2 . X„}  by  an  appropriate 

set  of  transformation  equations. 

The  easiest  transformation  to  implement  is  one  where  each  Y,  is  a  monotonic  function  of 
only  one  Xj,  i.e.. 


Y|  =  Yj(X|)  and  X,  =  X^Y>).  (3.59) 

Consider  the  integral  expression  (3.29)  for  the  two-parameter  model, 

T2(t)  =  H(xuxi,t)pXl(xl)pXl(xi)dxtdx1  (3.60) 

J  *2.1  J  *1.1 

which  represents  the  probability  of  survival  term  for  the  subspace  D2-  In  terms  of  the  new 
variables. 


where 


and 


Ut) 


'yia  r >1,2 

«//IJ'i,')Pt1()’iKi(),2M)’i<(»’2 

(3.61) 

y».\  J  j'i.i 

Py,(yd  -  Px,(xtyd)\dxJ<tytl 

(3.62) 

HJ.yiJi,0  =  ff(*i(>,i).*2(>’2).0 

(3.63) 

y*j  =  y{x,j). 

(3.64) 
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For  example,  if  y,  =  cptj,  where  c,  is  a  icaling  constant, 

H&iJi't)  =  HixJc^Jc^t),  (3.65) 

Pi, (yd  =  PxJiyt/ed/c,  (3.66) 

and 

yt.j  =  CiXu.  (3.67) 

A  particularly  useful  transformation  can  be  made  by  allowing  one  of  the  new  variables  to 
be  time  dependent  For  example,  consider  the  transformation  from  (X,.  X2, . . .,  XM}  to 
{Y„  Xlt . . .,  XM}  such  that 


Yt  =  Yj(X„l) 


X,  =  X,(Y,.f) 


and  the  transformation  functions  are  monotonic.  Again,  the  essential  mathematics  can  be 
illustrated  by  considering  the  Dt  probability  of  survival  term  for  a  two-parameter  model  satis¬ 
fying  the  assumptions  for  equation  (3.30).  The  transformed  expression  for  /,(x2,0,  (equation 
(3.31)),  is 


A(*a.O  = 


f>,.a( 

,0  = 

J  >...( 


where  p,,(y ,).  and  the  integration  limits  are  defined  by  equations  (3.62)  to  (3.64). 

The  expression  for  the  D2  contribution  to  Pr(t)  is  generated  by  taking  the  time  derivative 
of  equation  (3.30)  which  process  reduces  essentially  to  establishing  dlxjdt,  viz.. 


^  =  J  HAylsi,t)^Pil(yi)dyi+ 

+^jfpr,(yi 

(3  70) 

noting  that  the  second  term  in  (3.70)  results  from  the  fact  that  Py,(>i)  »  time-dependent. 

Before  proceeding  further,  it  is  necessary  to  examine  recasting  some  of  the  relationships 
implied  by  the  transformation  equations  such  as, 

ffrij'i.O  *  yi-yi(*i,‘)  ■  o.  (3.7i) 

Equation  (3.71)  can  be  interpreted  as  an  implicit  function  for  xt  in  terms  of  the  independent 
variables  y,  and  t.  Now, 

iF  iF  7>F 


which  lead  to, 


dx.  =  ^4y.+^dt 

*yr  W 


IlL  +  +  £*+  1^1*  -  0 

Wl  +  +  (W  +  W  J 

l1  -  if,  s^r1  ~ iir + irirr  ■ 


Because  y,  and  t  are  independent,  the  relationships. 


IS 


^i  =  l/^?  = 

fcXj  /  it  I  it 


follow. 


Returning  to  the  evaluation  of  dljdt, 


i  f  ,  .ThXiliXi  i  ixt  ._Jixt\ 

-  +'»'<*■>*  »  *nw 

a  r  /  »**i  i*x i  .  v**iYi 

-  “i“!(372)- 

The  second  term  in  (3.70)  can  be  integrated  by  parts 

fyt.2  d 

#/3’i.*a.<)r,/>v,0’i)<6'i 

J*,.  a< 

H£yirX1,t)pri(yiy^  I  +  I  — 7f>(.>’l,*2>fh^-PY,(J’l)<6’ 

J3'I.l  J  J'i.i  1 


so  that  (3.70)  becomes. 


dl.  Cy'*{i 

^‘L,r 

J  1H, 

l  *  11  ,J 

<6-1.1  toil  1„ 

— ar- arl  p 

r*, 

where  the  first  term  reduces  finally  to  I  — 

J  y«.i 


+  ^IH»(3'i’**>,)l"^r  Py,(>i)<6,i+ 


ItfO’i.jfi.OJMJ’iMvi- 


Equation  (3.73)  can  be  used  to  generate  the  complete  expression  for  p/t),  following  the  analysis 
leading  to  (3.35).  Examides  of  the  use  of  transformations  to  time-dependent  random  variables 
appear  in  the  following  sections  of  this  report. 


3i  Summary  of  Steps  Required  ta  Cm  irate  a  Reliability  Model  for  Fatigoe 

Having  detailed  the  derivation  of  the  reliability  functions  for  a  fatigue  model,  it  is  worth¬ 
while  summarising  the  essential  steps,  together  with  pointers  to  the  relevant  equations,  required 
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to  generate  a  reliability  model.  The  model  derivations  described  in  the  following  sections  of  the 
report  follow  the  steps  described  below. 

(i)  Random  variables  representing  parameters  in  the  fatigue  process  are  defined  and 
corresponding  density  functions  established. 

(ii)  The  various  time  zones  in  the  fatigue  process  are  identified. 

(iii)  For  each  time  zone  a  relationship  between  risk  rate  and  the  random  variables  is  estab¬ 
lished,  equations  (3.5)— (3.13). 

(iv)  Boundaries  between  the  time  zones  are  delineated  in  terms  of  the  random  variables, 
equations  (3.14),  (3.15). 

(v)  If  desired,  transformed  random  variables  are  defined  and  the  equations  established  to 
this  stage  recast  in  terms  of  those  variables. 

(vi)  An  integral  expression  for  P^t)  is  obtained  following  the  analysis  defined  by  equations 
(3.16H3.22). 

(vii)  The  integral  expression  for  P^t)  is  modified  to  indude  the  effects  of  inspections.  For 
the  models  considered  here,  the  modifications  amount  to  changes  in  the  positions  of 
some  of  the  time  zone  boundaries. 

(viii)  P/t)  is  defined  using  equation  (3.27). 

(ix)  />r(f)  and  r(t )  are  generated  by  taking  the  time  differential  of  P^t).  This  follows  the 
example  analysis  of  equations  (3.29>-(3.35),  using  (3.73)  if  transformed  variables  are 
used. 

(x)  The  integral  expressions  for  PJt)  and  p^t)  are  transformed  so  that  (3.38)  and  (3.40) 
can  be  used  to  evaluate  p^Rjr)  and  Pa(ff|F  >  t)  (or  corresponding  densities  for  some 
other  structural  characteristic). 

(xi)  The  effections  of  additional  risks  are  included  as  detailed  in  Sections  3.5.1  and  3.5.2. 


4.  THE  TWO-PARAMETER  MODEL  OF  PAYNE  ET  AL. 

One  of  the  functions  of  this  report  is  to  provide  a  theoretical  basis  for  the  documentation 
of  the  computer  program  called  NERF  which  evaluates  the  reliability  models  developed  by 
Payne  and  various  co-authors.  It  is  appropriate  that  the  analysis  of  the  previous  sections  be 
applied  first  to  the  derivation  of  one  of  these  models. 

The  reliability  model  derived  in  this  section  is  based  on  the  earliest  of  the  Payne  models 
and  relies  on  similar  assumptions.  The  only  significant  difference  is  the  treatment  of  the  virgin 
risk  for  uncracked  structures.  The  Payne1'4  analysis  derived  risk  rates  for  cracked  structures 
only;  the  risk  term  for  the  uncracked  time  zone  was  ignored.  This  term  is  included  in  the  present 
analysis  for  completeness  and  the  Payne  et  al.  expressions  can  be  obtained  by  setting  r,  =  0. 

The  expressions  developed  here  correlate  most  closely  with  those  presented  by  Payne  and 
Graham,4  although  reference  to  the  most  detailed  of  the  earlier  papers  (Diamond  and  Payne3) 
will  also  be  made. 

4.1  (mimiUnai  and  Model  Eqnodona 

41.1  Random  Variables 

The  model  assumes  that  the  variation  in  the  fatigue  process  can  be  represented  by  two 
random  variables. 

(i)  X, :  Comparative  fatigue  life 

Each  structure  in  the  population  is  assumed  to  have  a  fatigue  life  t,  given  by 

tf  ■  IfX ,  (4.1) 
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where  ft  is  a  known  median  fatigue  life.  A  suitable  density  function  />x,(*i)  is  also 
known. 

The  parameter  xt  is,  in  effect,  a  time  scaling  parameter  which  relates  the  time 
history  of  the  structure  to  a  median  time  scale.  For  example,  if  f,  is  the  median  time  to 
crack  initiation,  then  the  initiation  time  for  a  structure  with  fatigue  life  lt  is  xtf,. 

(ii)  Xj :  Relative  residual  strength 

For  any  structure  in  the  population,  the  ratio  of  the  strength  R  at  a  given  crack  length 
to  the  median  strength  R  at  the  same  crack  length  is  assumed  to  be  constant  for  that 
structure.  This  constant  is  the  second  random  variable  in  the  model 


= 


m 

R(o) 


(4.2) 


4.1.2  Time  Zones  and  Physical  Processes 

The  three  time  zones  indentified  in  Section  3.2  apply.  The  boundaries  between  the  time 
zones  are  given  by 


r,  =  xjt 

for  DJDi 

(4.3) 

h  =  Vr 

for  Z)j/Dj. 

(4.4) 

An  additional  condition  that  all  structures  have  a  strength  greater  than  a  certain  specified 
minimum  value,  Rmln,  also  applies.  Thus  structures  in  Dt  or  D2  must  satisfy, 

*  >  (4  5) 


4.1 3  Risk  Rate  Equations 

The  relationship  between  strength,  crack  length  and  time  is  assumed  to  be  known  in  the 
mean,  i.e., 

R  =  *W0)  (4.6) 


where  the  functions  have  the  general  form  shown  in  Figure  4.1.  For  a  structure  in  D ,  the  strength 
is  given  by 


H  —  X2Rq. 

(4.7) 

For  a  structure  in  D2  the  strength  is 

R  =  Xj^r/x,)). 

(4.8) 

Defining 

<K'/*i)  =  R(a(t/xi))lRo 

(4.9) 

as  the  median  relative  strength  decay  function. 

R  =  x2^0^(r/x,). 

(4.10) 

The  risk  rates  for  Z>,  and  P  are. 

r,  =  r,(x2£0)  m.  Pt(X2R 0)  .  1, 

(4.11) 

rx  =  ^(xj/Jo^O/x,))  e  PilXiRoWIxJ) .  1,. 

(4.12) 

The  subspace  boundaries  can  now  be  defined  in  terms  of  the  values  of  the  random  variables. 
For  structures  in  2>„ 
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x2  >  R-jIL  and  //f,  <  xt  <  oo. 


(4.13) 


For  structures  in  Z>2, 

xj  >  and  <  x,  <  r/f,  (4.14) 

where 

1U  =  ®in{/„*- *(J?_  J  U)}.  (4-15) 

or  in  the  reverse  order, 

lilt  <  x,  ^  r/f,  and  iW(A>W/Xi))  «  x,  <  oo.  (4.16) 

The  positions  of  these  boundaries  in  x  space  are  illustrated  schematically  in  Figure  4.2. 

4.2  TraMfonutkms  of  the  Random  Variables 

It  is  convenient  to  recast  the  model  in  terms  of  two  new  random  variables. 

(0  Vj:  Age 

V»  -  f/X,.  (4.17) 

Whereas  all  structures  with  a  given  value  of  X]  have  the  same  fatigue  life,  at  time  t 
structures  with  the  same  value  of  Y ,  have  the  same  value  of  relative  strength  decay.  In 
this  way  Yt  can  be  interpreted  as  a  measure  of  the  age  of  the  structure.  The  probability 
density  for  Y,  is, 

Pr,(yt) = PxtO/yMyl  (4-i») 

(ii)  Y2:  Virgin  strength 

Yj  =  Xa*o.  (4.19) 

All  structures  with  a  given  value  of  Y2  have  the  same  strength  before  crack  initiation. 
This  strength  is  called  here  the  virgin  strength.  The  probability  density  for  Y2  is, 

PvXHj)  =  P*0'j/A>)/A>.  (4-20) 

In  terms  of  the  transformed  random  variables,  the  risk  rate  equations  are, 

'i  -  ',(*)  (4  21) 

and 

'2  =  rt(y2tKyi)).  (4.22) 

For  structures  in  D„ 

A.i.  <  yi  and  y,  <  f,  (4.23) 

and  for  structures  in  Dlt  either 

jR.„  <  y2  and  ?,  <  yx  <  hs,  (4.24) 

where 

Its,  -  wn{l„rl(*m  Jyi)h  (4.25) 

or 

ltKyi<  l,  and  <  yx-  (4.26) 

The  positions  of  the  boundaries  in  j  space  are  illustrated  schematically  in  Figure  4.3. 

4J  An  Integral  Expresses  far 

Using  (3.2), 

P\IM  -  exp{-r,(x2/lo)t}, 


20 


(4.27) 
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Fj(x,f)  =  exp 


or  in  terms  of  y. 


Pfo,i)  =  expt-r,^)*} 


Pfaj)  =  exp 


X: 


I  'xiyi'Kt'yilW 

*i 


The  integration  with  respect  to  t'  can  be  transformed  by  letting 

y'i  =  t'yjt 

so  that 

dt'  =  Olyjdyi,  t'  =  t  «*>  >,'  =  y, 

and 

r<  ==  =  /,//*  «•  =  f,. 

The  expression  for  P«(y,f)  is, 


/»5(y,0  =  exp 


—  exp 


-ri(y2Wlyi 


r2(yi*(y\))dy[\ 


-( tlvi ) 


rstotf'O’OVO'i'  |- 


(4.28) 

(4.29) 

(4.30) 


(4-31) 


Note  that  the  term  in  square  brackets  in  (4.31)  is  independent  of  time.  This  fact  is  used  to  advan¬ 
tage  during  the  numerical  evaluation  of  the  reliability  functions  by  the  NERF  computer 
program  and  is  a  significant  motivation  for  the  transformation  from  X  to  Y. 

The  integral  expression  for  P^t)  can  now  be  constructed  using  equation  (3.20)  and  the 
limits  defined  by  (4.23H4.26), 


P&) 


f 

=  Pv,l 

J*... 


,0'a)exp{  -rl(y1)t}«!>'2i,vI(fi)+ 


+ 


f«>  f'r,. 

* 

\ 

1  Pv.fJ'j)  exp 

-(</>.) 

»'l(>’2)f|+ 

ri(yiHy{)W  j 

J  J  1. 

• 

f,  J) 

\dy2 


(4.32) 


The  early  models  developed  by  Payne  et  al.  assume  that  inspection  procedures  detect 
deficient  structures  according  to  a  crack  length  criterion.  The  inspections  are  assumed  to  be 
perfect  so  that  all  structures  with  cracks  greater  than  the  criterion,  denoted  hereby  a*,  arc  removed 
from  the  population.  In  terms  of  the  reliability  model,  this  means  that  at  the  moment  of  inspec¬ 
tion,  all  structures  that  have  aged  to  lt  where 

fs -«•*(«*)  (4.33) 

are  removed.  For  the  J\h  inspection  at  lij  say,  all  structures  with  xt  <  tijIU  are  removed  so  that 
the  effect  of  the  inspection  can  be  represented  in  equation  (3.21)  by  the  function, 

Xx,rtj)  -  HAxt-ttJfJ  (4.34) 

where  H£x)  is  the  unit  step  function  such  that  HJix)  —  0  for  x  <  0  and  H(x)  -» 1  for  x  >  0. 
The  corresponding  (Unction  for  yt  when  t  >  tt,  is. 
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(4.35) 


S(y,tij)  =  *.('/* 

All  structures  with  yx  >  lit/tij  have  been  removed  from  the  population. 

The  effect  of  an  inspection  at  tij  can  be  incorporated  in  (4.32)  by  adjusting  the  upper  limit 
of  integration  for  yv  This  is  achieved  by  replacing  I,  in  equation  (4.25)  by  if  where 

If  =  min{ff,fd///ij}.  (4.36) 

If  the  inspections  occur  very  frequently,  then  the  removal  of  structures  may  be  regarded  as 
being  continuous.  For  this  case, 

If  =  min{f„fd}.  (4.37) 

For  the  remainder  of  this  section  the  upper  limit  for  the  yx  integration  will  be  assumed  to  incor¬ 
porate  the  effect  of  inspections  prior  to  the  time  that  the  integral  expression  is  being  evaluated. 


4.5  Generation  of  p^t)  and  Risk  Rates 

Equation  (4.32)  can  now  be  differentiated  to  obtain  the  expression  for  Pf(r).  Writing 

W  =  TM+T^t) 

and  noting  that 

=  -i-  =  ^L 

it  X,  t' 


ML- 

-I 


[(y2>i(yi)eM-rt(yi)>}*iy2prt(i,)- 


Py,(yi)exp{-t‘i(y2)t}dy2 . fy,(W/, 

where  (3.73)  has  been  used  to  obtain  the  second  term.  Tz(t)  can  be  written  as. 


Tt(  t) 


r°°  r*< 

J  Rtmia  J  U 


HJ.yi,yi,t)pjl{yx)dyl 


where 


The  results. 


.0  =  cxpj 


j: 


-(*lyi)\ri(y1)li+  rz(y2<li(y't))dyl 


dHy^ 


dt 


Htft.n>y2>i) 


J)  =  expj— (r/ffjl)  ri(y2)f,+J^ 


riiyi+WWi 


and 


*  *xp{— r2(y2)f} 


may  be  used  to  yield, 


(4.38) 


(4.39) 


Ft/J'j) **  *«»»»  — (t/j'i)  '•|(y2)l,,+ ^(yW'WM'i' 


.X 


-J*^  ri(.y2)ti+ j*^  '*  ri(yi'Hy[)W\  Py, 

fco 

+  P\1(y&*Q{-ri(yi)t}(hlt)[>Yl(id<iyi-  (4.40) 

JXmi* 

As  expected,  the  second  term  for  dTJdt  cancels  with  a  similar  flux  term  for  dTJdt. 

The  second' term  in  (4.40)  represents  the  flux  of  structures  from  the  D2  subspace  and  can  be 
decomposed  in  the  following  manner.  Equation  (4.2S)  implies  the  following  relationships. 

ff.„  =  'l'~i(RmJyi)  =  U  say,  for  y2  <  f?) 

ffjJ  =  t*  for  <  oo. 

The  integration  over  y2  can  be  divided  into  two  sub-integrations.  The  first  follows  the  line 
/Li.  =  y2*Ky\)  (Fig.  4.3),  and  yields  the  total  flux  of  structures  that  fail  because  their  strength 
has  fallen  to  Rmin-  The  second  integration  yields  the  flux  of  structures  that  have  reached  their 
fatigue  life.  In  the  second  case,  if 

?*  =  Ittjtij  <  lt, 

(dy.  2  SVi  \ 

- —  I  in  (3.73)  is  zero  so  that  there  is  no  contribution  to  p^t) 

from  this  sub-integration.  This  corresponds,  of  course,  to  the  condition  that  a  previous  inspec¬ 
tion  has  removed  those  structures  that  would  have  reached  their  fatigue  life  limit  by  the  time 
that  Pr(t)  is  being  evaluated. 

Bearing  the  above  considerations  in  mind,  the  expression  for  p^t)  becomes. 


ri(y2>Kyi)}tX] 


p  -(t/yi) 


'•lO’iMH- 


r°° 

|  =  PYj(y2)r  i(3'i)cxp{ — r ,( yjV  HM>*|(0+ 

J  /L.. 

+J^  Pr^j^  '■i(yJlMyi))e*p|-Wy,)|r1(y2)fl+ 

1)  CRmlJ'Kt?)  , 

+  1  r2(y2i>(y{))dy'A\ihtt(yl)dyldy2+  I  PtM^xp -(///„)  r1(iy2)f,+ 

, 

+  J  ^  r2(y2\f>(y{))dy[  pr,(fK)dy1+ 

J  *»jm  '  { 


r2(yi'Hyi 


,(W »+ 


nwx 


Mtt  yfyi- 


(4.4iyt 


The  total  risk  (ri0=P»(0/F»(<)).  can  be  written  as  the  sum  of  three  risk  terms  as  in  equation 
(3.36).  These  risk  terms  are: 


Py,0'2>’i0'2)«p{  ~rt(yt)t}dy: 


■  I  Pv. 

Jo 


(y,Xy./F,(0. 


ti(2r^)-0  x+y, 
-1  x  =  y. 
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/*« 

r.(  0  =  Py  2(yi)  r2(y2\li(y1))H^yi,y2,t)pyXyi)dyi<fy2lp»(t) 

J  Pm\a  J 

rM=  p*JLyjjn/l»v»t)*tVMrMQ- 

J  ^min 


+#V?) 


r.r 

J  *.,./■ 


P\1(yi)1H,(it<yiJ)PYl{h)<iyilP*{‘)- 
\  • 


The  computer  program  NERF  evaluates  the  expressions  given  above  for  p^t)  and  the  three 
risk  rates.  For  direct  comparison  with  the  model  developed  by  Payne  et  a/.,1-4  the  risk  of 
failure  of  uncracked  structures  is  neglected  by  setting  rt  to  zero.  The  integral  expressions  for 
r,  and  rr,  in  terms  of  the  original  random  variables  xt  and  x2,  are, 


foo  ft/f, 

=  PxS.Xl) 

J  RmlJ^O  J  t fit. XI 

exp  -  |  r2(RoXi<Kt'i 


2(K0X2M<!Xl)) 


APxSx^idxj. 


Pxfcih — 
Xmd*0  UX1 


ff' 

l  J  *i« 


r2{R0x2Wff.x,lt))dt'\pxl(tu 


Equations  (4.4S)  and  (4.46)  compare  directly  with  equations  (10)  and  (1 1)  in  Payne  and  Graham.4 
The  only  difference  between  the  expressions  given  above  and  those  developed  by  Payne  and 
Graham  is  that  the  lower  limit  for  the  x2  integration  is  set  at  RmljK0  here  to  correspond  to  the 
known  limit  for  surviving  structures.  Payne  and  Graham  do  not  identify  the  two  sub-integrations 
for  the  rr  term  and  attribute  all  fatigue  fracture  failures  to  failures  with  R  —  RmiB.  Noting  that, 

lf  =  min  {f?.*  -  \RmJ Vs)},  (4.47) 

their  analysis  requires  that  lf  >  ^~i(RmijR0x2m„)  where  x2(MI  is  the  numerical  value  of  x2 
beyond  which  the  density  function  pXl{x 2)  is  insignificant. 

This  means  that  the  strength  function  (3.8)  must  be  known  for  strength  values  lower  than 
Hal,.  If,  as  is  often  the  case,  there  are  insufficient  experimental  data  to  define  the  strength  function 
for  very  low  values  of  strength,  then  it  may  be  preferable  to  truncate  the  function  and  impose 
a  limit  on  fatigue  life  via  ff.  The  two  r,  terms  in  (4.41)  will  then  exist  and  may  be  used  to  determine 
the  contribution  made  by  the  poorly  defined  section  of  the  input  data  to  the  total  risk  rate. 


4.6  Probability  Density  for  Strength 

The  integral  expression  for  Pg(t)  can  be  transformed  from  the  random  variables  Y1(  Y2 
to  Yt,  R  where 

Y2  =  R/lKYj)  (4.48) 


SQ'tO's)  1 
»(*,*)  "  Kyi) 
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so  that  (4.32)  becomes, 


r°° 

P*(t)  =  Pi1(R)cxp{-rl(R)t)dRPyt([l)+ 

J  Kum 

+J^  (4.49) 

Defining  the  strength,  R,  as  a  structural  characteristic  of  interest,  equation  (2.17)  can  be  used  to 
generate  the  conditional  density  of  strength  given  survival  to  time  t, 

/>«(rt|F  >  t)  =  [ /^(/tlexpf  ->■,(*)/  }^v,(A)  + 


V? 

PvJ 

J  *. 


,W»(J '.))-^l^y-)j>Y,(3'1V»i  / ^0.  (4.50) 


4.7  Failure  Density  for  Strength 

The  expression  (4.41)  for  p^f)  can  be  transformed,  using  (4.48),  into  the  form. 


MO  = 


foe 

+  f" 

J  *»,„  J  *. 

+J ^ 

+WtJ7)  PvJLm 

J  *«.. 


(R)exp{-rt(R)t}dRPri(!l)+ 


^*.>'2.0/>Y,(^)<fy2  + 


+W&  (4.51) 

where  the  third  term,  being  an  integration  along  the  line  Pid'O'i)  =  J?„|„  has  not  been  trans¬ 
formed.  The  resulting  expression  for  the  failure  density  of  strength  is, 

M*\t)  =  (»•,  WPr/^xpi  -r1(*)t}P¥l(f,)+ 


W«+ 


+**,*.J 


H- 


PvtyiYjH^WiyilOPvi  <>a>«  /  JWO 


The  last  term  (which  is  obtained  from  the  third  teim  in  equation  (4.51)),  represents  the  con¬ 
centrated  density  due  to  failures  with  Jf  =* 

Diamond  and  Payne1  derived  an  expression  for  the  probability  distribution  for  the  load  at 
failure.  Because  the  load  at  failure  is  equal  to  the  strength  at  failure,  this  distribution  is  equivalent 
to  the  cumulative  conditional  distribution  of  strength  given  failure  at  r. 
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=  f  M*\‘)dR.  (4.53) 

Their  expression  accounts  only  for  the  losses  by  static  failure  (r2)  and  makes  the  approximation 
(see  following  section),  that  the  loss  factor  Hy(yuy2,t)  «  1.  From  the  second  term  of  (4.52), 

r  r  rn 

/V(*|/)  *  I  r2(R)  I  '  Pr,(yi)p*^~~dytdRIMt), 

and  transforming  back  to  the  original  random  variables,  x„  x2,  with  the  additional  approxi¬ 
mation  that  Pa(t)  k  1  so  that  p*(t)  «  r(t), 

ft/?,  fR/iJtoWM) 

Pt(R |0  *  r2(RQx2il/(t/xl))pXl(x1)pTl(.X2)dJCldx1/r(t)  (4.54) 

J  W,  J  0 

which  is  the  same  as  the  expression  given  by  Diamond  and  Payne. 


4.8  Approximate  Expressions  for  the  Risk  Rate 


During  the  initial  development  of  the  reliability  models  numerical  difficulties  associated 
with  the  evaluation  of  the  integral  expressions  provided  a  strong  motivation  for  simplifying 
reliability  functions.  Accordingly,  approximate  expressions  were  derived  (e.g.  Refs.  1-3),  by 
recognising  that,  in  many  practical  examples,  interest  is  confined  to  the  early  part  of  the  popula¬ 
tion’s  time  history  where  Pa(t)  as  1  and  only  a  very  small  section  of  the  population  has  aged  to  its 
fatigue  life  limit  so  that  H  &  1.  Making  these  approximations  and  setting  /?m,n  =  0,  the  expres¬ 
sion  for  the  risk  of  static  fracture  by  fatigue,  r,(t),  becomes 

/•oo  n/r, 

r.( 0  *  Px2(*2)  r1(.KoX2\l>(tlx1))Pxt(Xi)dxldx1  (4.55) 

JO  J  tin 

which  corresponds  to  equation  (3.11)  of  Diamond  and  Payne.3  The  risk  of  fatigue  fracture  is 
given  by. 


'K0 


00 


0 


(4.56) 


Although  this  approximation  has  been  used  for  r^t),  an  alternative  was  derived  by  Payne.2 
Here,  this  approximate  expression  can  be  deduced  by  neglecting  the  effect  of  the  risk  in  D2. 
If  r2  (and  r,)  are  set  to  zero,  H(xux2,t)  =  1  and  the  expression  for  the  probability  of  survival  is. 


from  which. 


so  that, 


PdO  = 


dPdf) 

dt 


Too 

=  Pxt(xi)dxx 

Jt/I' 


-Pasum 


r(t) .  (=rKf))  =  (4-57) 

Px  i(*l)<&l 

J  */!, 

This  expression  is  in  agreement  with  that  derived  by  Diamond  and  Payne,3  with  the  exception  of 
the  factor  l/ff  which  was  omitted  in  Reference  3. 
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5.  THE  TWO-PARAMETER  MODEL  OF  HOOKE 

The  reliability  model  developed  by  Hooke5  “*  uses  two  random  variables  and  is  closely 
related  to  the  Payne  model  described  in  previous  sections  of  this  report.  In  fact,  Hooke’s  model 
can  be  considered  to  be  encompassed  by  equation  (4.32)  for  P^t)  which  is  more  general  than  the 
corresponding  expressions  derived  by  Payne  or  Hooke.  However,  rather  than  simply  enunciate 
the  Hooke  equations  directly  from  equation  (4.31),  Hooke’s  model  will  be  derived  here  directly 
using  the  approach  summarised  in  Section  3.6. 


5.1  Assumption!  and  Model  Equations 

Xt :  Characteristic  time. — Hooke  lets  an  appropriate  event  define  a  characteristic  time  for 
a  structure,  and  makes  the  assumption  that  all  events  in  the  history  of  a  structure  are  related  to  a 
median  time  history  by  a  single  definition  of  characteristic  time.  This  means  that  xt  is  the  same 
time  scaling  parameter  defined  by  equation  (4.1)  for  the  Payne  model. 

Y2:  Virgin  strength.— The  second  random  variable  in  the  Hooke  model  is  the  same  as  the 
virgin  strength  parameter  defined  by  equation  (4.19),  following  the  definition  of  relative  residual 
strength,  for  the  Payne  model. 

(Note  that  although  the  virgin  strength  is  a  basic,  as  opposed  to  a  transformed,  random 
variable  for  the  Hook  model,  it  is  denoted  here  by  Y2  to  emphasize  the  equivalence  with  the  trans¬ 
formed  variable,  Y2,  in  the  Payne  model.  As  far  as  is  practicable,  the  notation  developed  in  Section 
4  will  be  used  in  this  section.) 

The  model  developed  by  Hooke  corresponds  to  one  in  which  the  third  time  zone  does  not 
exist.  There  is  no  imposed  fatigue  life  limit  and  no  specified  value  for  minimum  strength.  The 
boundary  between  Dx  and  D2  is  given  by. 


<1  =  Xj{. 

(5.1) 

The  model  equations  are  similar  to  those  described  in  Section  4.1.3,  so  that 

<T  =  't(j Ks) 35  Piiyz)  ■  4 

(5-2) 

and 

r2  =  r2(y2^(f/x,))  s  A.(>i*(»/*i))  •  4- 

(5.3) 

For  structures  in  Z>„ 

0  <  y2  <  oo  and  »/f,  <  x,  <  oo 

(5.4) 

and  for  those  in  D2 

0  <  y2  <  oo  and  0  <  xt  <  t/f,. 

(5.5) 

53  An  integral  Expression  for  PJf) 

Using  the  model  equations, 

Pt(xx,y2,t)  =  exp{-r,(yj)r} 


and 


PfaiJj't)  -  exp  -r,(y»)x,f,-  J  rj(yX»7*i)>*'  • 

J  xJt 


The  expression  for  P^t)  is, 

iVM- 


J  0  J  III, 


e*Pi-ri(yi)t}Pxt(xl)r*,t(y1)dxtify2+ 


(5.6) 


(5.7) 
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Writing, 


■J.T 

■J: 


e*P 


ri{yiW!xJ)dt'\pxXxl)l>v1{y2)dxldy2. 


>1(^2).  t  <  xji 

'XyiHl/Xi))  = 

ri{ytH‘lx  1)),  t  >  x ,f„ 

the  two  terms  in  (5.9)  can  be  combined,  to  produce 


(5.8) 


(5.9) 


*vo= 


r'(y2Wi'lx  ,))<*' 


Pxl(xl)pYt(y2)dxidyi 


(5.10) 


which  is  identical  to  the  expression  obtained  by  Hooke  (equation  17  in  Ref.  8).  Note  that  the 
combined  expression  does  not  explicitly  identify  the  contribution  made  by  the  failures  of  un¬ 
cracked  structures  but  relies  on  that  contribution  via  the  definition  of  r'. 


5.3  Generation  of  Pr(t)  and  Risk  Rates 

Equation  (5.10)  can  be  differentiated  to  obtain  the  expression  for  p^t). 


’00  r ao  (  ft 

MO  =  '■'(>’2<Kf/*i))expj  -  rXy^it'IXiW 

J0J0  (  Jo 


Px,(xi  )Pn1(y2)dx1dy2  (5.11) 


The  resulting  expression  for  total  risk  (M0//g(0),  is  identical  to  equation  (19)  obtained  by 
Hooke.8 

Alternatively,  the  total  risk  can  be  considered  to  be  the  sum  of  two  risks  representing 
failures  in  uncracked  and  cracked  structures.  Differentiating  (5.8)  yields. 


and 


where 


f  00  Too 

r»(0  —  I  ,i(3,2)cxp{  -ri(y2)t}p*l(y2)dy2 1  />*,(*, V^i/MO 

Jo  J  f//, 

'2(y2'l'(t/xt))H(xl,y1,t)pXl(xl)pv2(y2)dxldy2/Pt{t) 

0 


H(xlty3,t)  =  expj 


-'’i(>2>*ifi-  f  r2(.y2'Kl'IXi)dt'\. 

J  Vt  j 


(5.12) 


(5.13) 


(5.14) 


5.4  The  Effect  of  Fatigue  Life  Limiting 

The  analysis  leading  to  equation  (5.10)  and  (5.1 1)  assumes  that,  there  is  no  imposed  limit  on 
fatigue  life.  As  was  done  in  section  4,  the  effect  of  such  a  limit  can  be  accounted  for  by  defining 
the  boundary  between  the  “cracked”  and  “failed”  time  zones  by, 

h  =  Mr-  (5.15) 

For  a  structure  with  X,  =  xu  equation  (5.15)  determines  the  fatigue  life  limit,  or  “point  of 
instantaneous  strength  decay”  as  identified  by  Hooke.  The  boundaries  of  the  D2  subspace  ate 
now  defined  by, 
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so  that  from  (S.10) 


0  <  >2  <  oo  and  t/i,  <  xt  <  t/l. 


(5.16) 


PM)  =  I  cxpj—  j  r\y%'Kt'lxi)W 
JO  J  t/l,  {  Jo 


Px,(xl)py1(y1)dxldy1 


(5.17) 


and  differentiating  (5.17) 


=  ''0,2lW*/*i))«p(- |  r'iytWIXiW 

JO  Jl/l,  1  Jo 

+“P  JQ  P»1(>a)cxp|-JVo'2^70)*'j 


Pxl(xl)pyl(yi)dxldy1 


JVC 0 


+ 


I  dy2 

PM)' 


(5.18) 


The  second  term  in  (5.18)  is  the  risk  of  fatigue  life  exhaustion.  The  corresponding  term  obtained 
by  Hooke  resembles  that  in  (5.18)  but  contains  errors  in  the  exponential  factor.  The  correct 
expression  in  his  notation  is,* 


ftCO  =  ferKO 


m(U0at.k,IHt)ltr0)dt 


f(U0)dU0IR(t). 


(5.19) 


Hooke’s  consideration  of  fatigue  life  limiting  concerns  only  the  r£t)  term  and  the  corresponding 
expressions  for  PM)  and  the  first  term  in  (5.18)  are  not  given.  If  fatigue  life  limiting  is  applied 
and  the  risk  of  fatigue  life  exhaustion  is  calculated  separately,  the  integration  limits  in  the  expres¬ 
sions  for  PM)  and  r(t)  must  explicitly  reflect  the  existence  of  the  boundary  between  the  sub¬ 
spaces  of  cracked  and  failed  structures  to  prevent  double  accumulation  of  failures. 

Returning  to  the  more  general  model  developed  in  Section  4  as  a  basis  for  the  discussion 
of  the  Payne  model,  the  />(/)  term  obtained  here  is  equivalent  to  the  second  term  in  (4.44)  with 
Rmi.  —  0.  In  practice,  the  motivation  for  applying  a  limit  on  fatigue  life  is  often  the  result  of  an 
inability  to  adequately  define  the  strength  function  for  low  values  of  ft  and  l,  will  correspond 
with  a  minimum  value  of  strength.  It  is  preferable,  therefore,  to  use  the  more  general  model 
which  can  be  reduced  to  the  Payne  or  Hooke  versions  by  setting  ?f  =  oo  or  Rmt,  =  0 
respectively. 


5.5  Failures  Caused  by  Fatigue  Weakening 

Consider  a  population  of  structures  which  do  not  weaken  by  fatigue.  The  only  non-empty 
subspace  is  D,  and 

r'  =  r,  =  rt(y2)  =  A.(Fi)  •  4-  (5.20) 

The  probability  of  survival  is  given  by. 


PM) 


-n 


00 


e*P{  -ri(yt)t  }Pxl(xl)pri(y2)dxtdy2 


=  I  *M~ri(yi)t}Pi,(yi)dy2 

0 


(5.21) 


which  is  identified  here,  in  line  with  Hooke,  as  the  survivorship  with  strength  preserved.  It 
follows  that  the  total  risk  which  will  be  called  the  risk  of  ultimate  failure,  ra(r),  is  given  by. 


rM'I? 


r  iO'i)*xp{  -r1(y2)f  JPY.O'aV&'j/JKO 


(5.22) 


*  This  notation  is  not  defined  in  the  table  of  notation  as  it  is  local  to  this  equation  (5.19)  only. 
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Returning  to  the  population  in  which  fatigue  weakening  is  occuring,  Hooke  identifies 
failures  that  have  occurred  because  of  the  fatigue  process.  Such  failures  are  those  of  cracked 
structures  by  loads  that  are  less  than  the  initial  (uncracked)  strengths.  Hooke  then  computes 
the  “probability  of  failure  of  weakened  structures  by  loads  below  their  virgin  strengths”,  Pf(t) 
say,  where 

PwO)  =  P&D-PJLt)  (5.23 

PS(t)  is  the  probability  of  survival  defined  by  equation  (5.21)  and  Pt(t)  is  defined  by  (5.10)  (or 
(5.17)  if  fatigue  life  is  limited). 

Hooke’s  approach  is,  in  fact,  incorrect.  Using  (2.6),  equation  (5.23)  becomes, 

PRO-MO-PXO  (5.24) 

which  on  differentiation  with  respect  to  time  yields, 

M‘)  =  jf(t)+PX  0-  (5.25) 

While  it  is  true  that  p/t)  can  be  decomposed  into  two  additive  terms,  that  term  representing 
failures  by  loads  greater  than  the  virgin  strengths  is  not  equal  to  pi(t)  (as  obtained  by  multiplying 

(5.22)  by  P&(t)).  Any  of  the  reliability  functions  derived  from  equation  (5.21)  cannot  account 
for  the  fact  that  failures  contributing  to  />?(/)  reduce  the  effective  population  of  structures  that 
can  fail  by  loads  greater  than  the  virgin  strengths.  pjj(<)  will  therefore  be  an  over-estimate  for  the 
contribution  to  p^t)  from  failures  that  have  not  occurred  because  of  the  fatigue  process,  and 

(5.23)  will  over-estimate  PJ(t). 

The  correct  approach  is  to  introduce  the  distinction  between  the  two  types  of  failures 
during  the  definition  of  the  model  equations.  This  can  be  done  for  the  Hooke  model  by  replacing 
equation  (5.3)  by, 

r2  =  'l(y2<l'(t/xi))+rt(y1)  (5.26) 

where 

rSUuKtM)  =  PitoHK*lxd)-PdyM  ■  4-  (5.27) 

For  structures  with  X,  =  xt  and  V2  =  y2,  r2  is  the  risk  rate  for  failures  by  loads  less  than  the 
virgin  strength,  y2.  The  expression  for  the  probability  of  survival  (with  fatigue  life  limiting)  is. 


■a 


exP{ -'■l(y2)t}pxXx1)pyJ(y2)dxldy2 + 


n‘lh 

H(xl,y2,t)pxXxl)pYl(y2)<ixi<ty2 
' /IS r 


where  the  loss  factor  in  D2  is  now. 


H(x2 j-j.t)  =  exp  -r2(y2 


=  exp{ 


-ri(y2)t-  r{{ 

J  xj, 

~r  i(  J's)*  }exp| — 


nnym’Mw 


r”(y2Wlx  dw 


The  expression  for  pjt)  is. 


'-n 


ri(y1)txj>{-rl(y2)t}px,(xi^tl(yi)dxldy2+ 
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Jo  J  t/i 

PXiOflt)  f  00  « 

-  “7 —  Pit 

Jo 


lti(y2)+ri(y2'lf(f/xl))]H(xlj2,t)pXl(xl)dxldy2+ 


t(y2)H(tl!t,y2,t)dy2 
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which  can  be  rearranged  to  yield, 

M0=  Pxt(xt)dxi  + 

JO  LJ  tlh 


f*  J-P  * 

J  t/l,  \  J  liXi 
r  oo  rt/i, 

PvI(yi)  r$ 

Jo  J  til, 

— I —  Pjl 

Jo 


,  Pxt(x  i)dx. 


) 


dy2+ 


r2(yi'Ktlxi))mxlj2,t)Pxl(xi)dxl<ly2+ 


{yi)H(tll,^2,t)dy2. 


(5.31) 


The  last  two  terms  in  (5.3 1)  represent  the  required  density,  pS(t),  for  failures  due  to  fatigue  weaken¬ 
ing.  The  first  term  is  the  density  for  failures  by  loads  greater  than  the  virgin  strength.  The  correct 
version  of  equation  (5.25)  for  this  model  is, 


i: 


Pr(t)  =  Pr(t)-Pr(t)+  />Y,(j’j)riO's)e*  p{  -'ite)*}  1  - 


m 

P*. 

Jo 

J  |l-expj-J  rUyWIxW  jpx,(*i)<ki 


(xjjdx,-!- 


(5.32) 


6.  TWO-PARAMETER  MODEL  OF  FORD 

Whereas  Payne  and  Hooke  derived  their  reliability  models  using  techniques  which  bear  a  dose 
resemblance  to  those  used  here.  Ford  solved  a  continuity  equation  for  the  probability  density 
for  crack  length  to  generate  reliability  functions.  The  functions  for  the  two-parameter  model 
for  “coarsely  random  cracking”  are  derived  here  and  provide  verification  that  the  two  different 
methods  yield  the  same  results. 


(.1  Assumptions  and  Model  Equations 

The  model  is  based  on  two  random  variables. 

(i)  Xj :  Comparative  crack  growth  rate 

Ford  assumes  that  the  equation  for  crack  growth  can  be  written  in  the  form, 

da/dt  —  X|  /(a)  (6.1) 

with  X|  =  1  being  equivalent  to  the  known  median  equation.  Equation  (6.1)  can  be 
transformed, 

da/dO/x,)  =/(o),  (6.2) 

so  that  x,  can  be  interpreted  as  a  scale  factor  for  time.  An  event  occurring  at  l—l0, 
where  f0  is  some  suitable  reference  time  in  the  median  time  history,  will  occur  at 
(l—lo)lxt  for  a  structure  with  X]  =  x,.  Clearly  X,  is  the  reciprocal  of  the  fatigue  life 
parameter  defined  in  Section  4  for  the  Payne  model. 

(ii)  X2 :  Crack  initiation  time 

An  obvious  choice  for  a  reference  event  is  crack  initiation.  If  this  event  is  assigned 
variation  and  a  random  variable,  X2,  an  event  occurring  at  time  f  in  the  median  time 
history  will  occur  at  time  t  in  the  history  of  a  structure  with  Xt  =  x,  and  X2  =  x2 
where, 

/ »  xa+(f — fi)/x«.  (6.3) 
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x2 


FIG.  6.1.  POSITIONS  OF  SUBSPACE  BOUNDARIES  IN  (X,.Xt)  SPACE 
FOR  THE  2  PARAMETER  FORD  MODEL. 
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Ford  imposes  fatigue  life  limiting  via  a  runaway  crack  hypothesis  and  because  there  is  no 
variation  in  strength,  a  limit  on  fatigue  life  defines  a  corresponding  minimum  value  of  strength. 
The  boundaries  between  the  three  time  zones  are, 

/i  =  x2  for  Dt/D2  (6.4) 

and 

h  =  (ff-fi)/*i+*2 

=  ft lx2+x2,  say.  (6.5) 

Since  there  is  no  variation  in  strength,  the  risk  rate  in  Dt  is  constant  and  was  included  by 
Ford  with  an  independent  “hijack”  risk  in  an  additional,  constant  risk  term.  In  this  analysis  the 
hijack  risk  is  not  included  during  the  development  of  the  model  but  can  be  included  subsequently 
using  the  general  results  obtained  in  Section  3.5.1.  The  model  equations  are, 

r»  =  Pl(& o)  •  K  (6.6) 

and 

r2  =  r2(xj(/— x2))  =  PM(lt  +*,(/  -x2))).  (6.7) 

For  structures  in  />,, 

l  <  x2  <  oo,  0  <  Xj  <  oo  (6.8) 

and  for  those  in  D2 

0  <  x2  <  /,  0  <  Xj  <  fi/(l-x2).  (6.9) 

The  locations  of  the  three  subspaces  in  (x2,x,)  space  are  shown  in  Figure  6.1. 


U  An  Integral  Expression  for  PJt) 

Using  the  model  equations, 

=  exp{— r,t}, 


/‘•(x.O  =  expj-riXj-J  r2{xx{t'-x2))dt'  , 


and  the  expression  for  P^t)  is, 


(6.10) 

(6-11) 


noo 

exp{-r,/}/»Xl(x1)p,J(x2>ix1<fc2  + 

r»  cmt-xj)  ft  ) 

+J  |  ^  exp  — rtjca— J  ri(xi(t—x2))dt']pxi(xl)px1(,xt)dxldxi 

-  exp{-r,t}^*l(t)+ 


ft  fti/(t-x2)  (  ft 

+  exp{-r1x2}/>*1(x2)  I  exp  -  r2(xx(t' -x2))dt'  PxS.xx)dx2<ix2. 

JO  JO  l  Jx2 

(6.12) 


43  Catenation  efpyft)  and  Blsfc  Rates 

Equation  (6.12)  can  be  differentiated  to  generate  the  expression  for  the  density  of  the  time 
to  failure. 


35 


Prit)  =  r1exp{-r,t}/'Xl(t)+ 


+  exp{-r1Jc2}pXl(x2); 

Jo 

=  rlexp{~rit}PXl(t)+ 


’  fhl0-x2) 

[‘  ) 

exp  - 

'i(*i(»' -x2)W iPxAxiidXi  • 

[Jo 

J  x2  j 

[t—x2)dx2. 

hKt-x  2) 


riMt-Xitexpl  - 


-x2))dt’\Pxi(Xi)dx  1+ 


Chl(t-x2)  (  ft  ) 

<Ht-x2)=  I  r2(xi(f-x2))expj-  I  r^x^t' -x2))dt’]PxS.xt)dxi  + 

,  ft 

+  77  1  TiexP  -  rAW  -x2)l(t-x2))dt'  pXl(?if(t-x2)) 

{t~Xl)  Jx2 

which  by  changing  the  variable  of  integration  in  the  exponential  term  becomes, 

f?i/('-*2)  |  fxi(t—x2)  ) 

<Kt—x2)  =  r2(Xi(f-x2))exp^  -  —  r1(c)dc)pXi(x1)dxl+ 

Jo  [  Xl Jo  J 


fhl(t-x2) 

J  fXi(t—X2)  ) 

r2{xi(t—x2))cxp  - 

—  r2(c)dc  Jj 

Jo 

XlJ  o  j 

(t-x2f 


exp  - 


(*-x2) 


ft i 

r2( 

Jo 


(c)dc}pXi(hl(t-x2)). 


Equations  (6.13)  and  (6.15)  compare  directly  with  (4.3)  and  (4.4)  presented  by  Ford." 
Apart  from  a  small  error  in  Ford’s  second  term  of  (4.4)  arising  from  an  algebraic  error  in  his 
equation  (3.4),  the  expressions  obtained  via  the  two  different  methods  are  identical. 

The  expressions  for  the  three  components  of  the  risk  rate  are, 

r,(t)  =  rlexp{~rit}PXi(t)IPn(t),  (6.16) 

C‘  Chl(t-x2) 


r,0)  =  J^exp{-rlx2}pXj(x2)J 

(  1  f  *i('-**) 


/*r 

i  x> 

xexp  -  - 

x\  o 


I  0 

r2(c)dc 


r2(xt(t-x2))x 

pXl(xl)dxldx2 

W) 


I’M  =  exp{  -r,x2}/)Xj(x2)~-  — ,  x 

n 


•m: 


r2(c)dc 


PxSf^t~x2))dx  2 


6.4  Trurfonwd  Random  Variables 

Having  established  the  equivalence  between  the  model  described  by  equation  (6.12)  and 
that  derived  by  Ford,  it  is  convenient  for  the  purposes  of  further  discussion  and  development 
to  introduce  new  random  variables  which  are  transformations  of  the  original  variables  used  by 
Ford.  Specifically,  the  introduction  of  a  random  variable  representing  age  will  facilitate  com¬ 
parisons  with  the  models  of  Payne  and  Hooke. 

Accordingly,  let 

Y2  =  X,  (6.19) 
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and 


Xj  >  t 

Y,  =  (6.20) 

/.+Y,(/-X2),  X2  <  t 

define  a  transformation  from  X  to  Y  with  Yt  being  the  required  age  variable.  For  X2  <  /  (or 
Yj  >  /,),  the  random  variables  Yt  and  Y2  are  not  independent;  the  joint  density  function 
Fv, . void's)  »  given  by, 

Pni.y1(yij2)  =  />i,(>2)/>x2('-(>t -h)lyi)lyi 

s  Pv1(yi)PvS.yi\yi)  (6.21) 

where  PY,(yt|y2)is  the  conditional  density  for  Y,  given  Y2  =  y2.  For  X2  <  t,Fyi(>,),  the  marginal 
density  for  Yt  is  given  by, 

r°° 

Prl(y\)=  PY1(y2)PYl(yi\yi)dy2-  (6.22) 

Jo 


For  Xt  >  /, 

PnXyi)  =  PxXWyiW/yl  (6-23) 

In  terms  of  the  transformed  random  variables,  the  model  equations  are  given  by  (6.6)  for 
r2  and 

'2  =  r2(y ,)  =  PM(y,))-  (624) 

(Note  the  difference  between  equations  (6.7)  and  (6.24)  in  the  definition  of  the  argument  for  r2. 
In  (6.24)  the  argument  of  r2  is  yx ;  in  (6.7)  the  argument  is  y,  —  f,.) 

For  structures  in  Dj, 

0  <  y  <  f„  (6.25) 

and  for  those  in  D2 

0  <  y2  <  00,  /,<>!<  min{fr,/,+y2i},  (6.26) 

as  illustrated  in  Figure  6.2.  Note  that  the  Ford  model  does  not  include  structures  with  x,  <  0 
(i.e.  structures  with  pre-existing  cracks).  This  condition  restricts  the  maximum  value  of  yt  for  a 
given  crack  growth  rate,  y2. 

Using  equations  (6.6)  and  (6.24), 

Fi(y,0  =  exp{-r,f},  (6.27) 


t-{yi-l\)!yi 


p»( y,0  =  exp|  -(/-(>!  -ii)lyt>t—  r 

J  t-iy^-ldbi 

1  f* 

(<-(>1  -fdlyiPi - r2(yl)dyl 

*J  h 


rz(yt+yz{t'-t))dt' 


Noting  that. 


nii  r » rf, 

exp{-rtt}prt,rt(yiJ'z)dyidy2  =  expi-rtl}  \  Ft1(>2)Pt,(J'iW>'i‘0'2 

0  Jo  Jo 


e*p{-r,»}FTl(/t), 


the  expression  for  the  probability  of  survival  is, 
W  =  exp{-r,<}FT|(fJ-(- 


+ Jo  Frito) **  Pr,(yt  ~ldlyz>i  ~ 

fz(yi)dy^iyidyz 
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FIG.  6.2.  POSITIONS  OF  SUBSPACE  BOUNDARIES  FOR  THE  EXTENDED  FORD  MODEL. 
SHADED  AREAS  DENOTE  INITIALLY  CRACKED  STRUCTURES 

(a)  (Xj,  X,)  space. 

(b)  (Y„  X,)  space,  (X,  —  initial  age) 

(c)  (Y„  Ya)  space. 


where, 

ltj2  = 

The  density  function  for  the  time  to  failure  is  given  by, 

MO  =  rie\p{-r1t}Plll(ld+ 


Myl)dy ,'WMy2  + 


f°°  , 

+  PvM)  »-i(j'i)Pv1(>,i|>,2)exp  -r,  .(/-(>» -f,)/y2)- 

J  o  Jr, 

~j2Jf  MyMy}^yxdy2+ 

+  f  yiPv&ilPySJbi'&V  I  My',)dy0dy2  (6.31) 

J  hit  yiJ  r,  ) 

which  can  be  readily  verified  to  be  equivalent  (but  based  on  a  reversed  order  of  integration),  to 
the  expression  given  in  equation  (6.13). 


6.5  Ford’s  One-Crack  Model 

Ford9  initially  derived  a  single  parameter  model  using  initiation  time  as  the  random 
variable.  This  model,  referred  to  here  as  Ford’s  “one-crack”  model  (although  strictly  all  models 
described  herein  are  for  a  single  crack)  can  be  deduced  from  the  equations  obtained  in  Section 
6.4  by  setting  Y2  =  1,  i.e., 

P*,(yi)  s  Py,(>'i|Y2  =  1) 

=  Px1('+f|->’i)>  (6-32> 

/>,(/)  =  exp{— rl/}/>Yl(/,)+ 

rmin{f„f,+/}  fn 

+  P¥,0'i)exp  Or,-  My'iWi  dy,  (6.33) 

Jr,  Jr, 


MO  =  '■1exp{-iy}/’Y,(f,)+ 

/*min{ff,r,+/}  (  Cy,  \ 

+  1  ^iK.^i^xpj -('+?, -I  ^  rdyQdyija 

+Hlt,t+t j!Pv,(ff )exp  — (/— fi)r1— r2(y,')dy1j. 


Equation  (6.34)  is  equivalent  to  equation  (4.7)  of  Ford. 


6.6  InspectkMs 

Ford10  presents  a  general  analysis  for  the  effects  of  arbitrary  inspections  on  his  one-crack 
model.  Here,  the  effects  of  perfect  inspections,  analogous  to  those  modelled  in  Section  4.4,  are 
included  in  the  two-parameter  model. 

The  effect  of  an  inspection,  at  tij  say,  is  the  removal  of  all  structures  with, 

xi{lij-x2)  >  rt-l,  (6.35) 

where  l4  is  the  life  corresponding  to  the  inspection  criterion  (equation  (4.33)).  At  time  t  >  tiJt 
the  effect  of  the  inspection  can  be  represented  by  the  removal  function, 

=  HJJt -r,+Xl(x2 - tij ))  (6.36) 
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or 


S(y,tij)  =  HXh-yi  +yi(t-tij)).  (6.37) 

These  removal  functions  can  be  either  inserted  in  the  integrands  of  the  reliability  functions 
or  represented  by  modified  integration  limits.  In  terms  of  the  transformed  random  variables, 
the  effect  of  the  inspection  at  tij  can  be  included  by  replacing  lt  in  (6.30)  by  If  where, 

It  =  min{ff,/d+j'2(/— tij)}.  (6.38) 

The  third  term  in  (6.31)  will  exist  only  when  If  =  I,. 


6.7  Failure  Density  for  Crack  Length 

Ford9  obtains  the  failure  density  for  crack  length  for  the  one-crack  model.  Using  the 
analysis  of  Section  2.6,  this  density  can  be  obtained  by  transforming  the  integral  term  in  (6.34), 

MO  =  r i(exp{  —rxt}P ¥l(a‘  '(ai))+ 


+ 


Of 


r2(a~  l(a))pYl 


(fl~‘(a)) 


+<Kff,t+f,)p»,(a',(af))cxp|  -(t-h)rl - 

where  a  =  a(t),  a,  =  a(/,),  a,  =  u(ff)  and  a'  =  dajdt. 
For  cracked  structures,  with  a,  <  a  <  a,, 


(6.39) 


pMO  =  -~r—Prt(o'  '(a))exp\ 


-0+1, -a 


(6.40) 


which  is  equivalent  to  equation  (4.12)  presented  by  Ford.' 


6.8  Model  Extensions  to  Include  Initial  Cracking 


As  postulated.  Ford’s  assumptions  do  not  allow  for  the  possibility  that  structures  may  be 
cracked  at  initial  time  (/  =  0).  Diamond  and  Payne3  presented  analysis  for  a  population  in  which 
all  structures  are  initially  cracked.  In  this  section,  Ford’s  model  is  extended  to  account  for  such 
structures  which  may  form  only  part  of  the  total  number  of  structures.  The  Diamond  and  Payne 
model  emerges  as  a  limiting  case  of  this  extended  model. 

A  structure  that  is  cracked  when  t  =  0,  can  be  considered  as  having  a  negative  crack  initiation 
time  and  commences  life  at  t  =  0  with  an  “initial  age’’  given  by, 

y,(0)  =  h-xlXl.  (6.41) 

Such  structures  can  be  included  in  the  Ford  model  by  extending  the  range  of  the  random  variable 
X2  to  (—00,00)  so  that  (6.26)  is  replaced  by, 

0  <  y2  <  co,  l,  ^  y,  <  l„  (6.42) 


(6.28)  by, 


P»(y,0  =  exp 


— max{0,f  — (y,  -f,)/y2)}ri  - 


_lf* 

*Jma 


My[W, 


and  (6.29)  by 


(6.43) 
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*s(0  =  exp{— r1/}i>Y,(^i)+ 


fco 

+  PY& i) 

Jo  J 


/>V,(>2bl)eXP 


— max{0,f— (.y,  -fi)/yj)}r,  - 


y% 

max{fi>y1-y2f} 


r2(yl)dy'i 


(6.44) 


Corresponding  equations  for  the  derived  reliability  functions  follow  trivially  from  equation 
(6.44). 


6.8.1  Initial  Age  as  a  Random  Variable 

The  model  based  on  equation  (6.44)  requires  that  initially-cracked  structures  are  represented 
by  deduced  negative  initiation  times.  A  more  convenient  random  variable  for  such  structures  is 
the  initial  age  defined  by  equation  (6.41).  Using  X3  to  denote  the  age  of  a  structure  at  /  =  0,  a 
structure  which  is  cracked  at  /  =  0,  has  age  Y,  at  time  t  given  by, 

Yl  =  X3+X,/  (6.45) 

where 

X3  =  f,-X2X2.  (6.46) 

The  definition  of  initial  age  for  an  uncracked  structure  is  somewhat  arbitrary.  One  definition 
which  allows  a  one-to-one  mapping  between  initial  age  and  crack  initiation  time  is, 

Yt  =  X,-H  (6.47) 

where 

X3  =  f,-X2.  (6.48) 

This  definition  amounts  to  the  assumption  that  the  variation  in  crack  initiation  time  arises  from 
variations  in  pre-aging  of  the  structure  rather  than  variations  in  the  damage  rate  (as  implied  by 
equation  (6.20),  which  assumes  all  uncracked  structures  have  zero  age  at  t  —  0). 

The  reliability  model  can  be  expressed  in  terms  of  the  random  variables,  age  (Y,)  and 
comparative  crack  growth  rate  (Y2), 


Y2  =  X„ 

(X 3+/,  Y,  <  /, 

(6.49) 

Y>  = 

(x3+y2/,  y,  >  i , 

(6.50) 

so  that  for  structures  in  Du 

Pyt(yi)  =  P\>(y\-i) 

(6.51) 

and  for  those  in  Dz 

PiyiyMi)  =  Px^yi  -yjO- 

(6.52) 

Equation  (6.44)  can  then  be  used  to  determine  the  probability  of  survival. 


6.8.2  The  Initial  Crack  Model  of  Diamond  aad  Payne 

The  initial  crack  model  developed  by  Diamond  and  Payne  can  be  deduced  from  the 
equations  developed  in  the  previous  section.  Assuming  that  all  structures  are  initially  cracked,  it 
follows  that  Dt  is  empty  and  that  D%  is  defined  by, 

0  <  y2  <  hji,  h+yit  <  yt  <  h  (6  53) 

(see  Fig.  6.2). 
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The  equation  for  the  probability  of  survival,  (6.44),  becomes,  using  (6.52), 


Pdf)  = 


'h H 

Jo 


PxXyi) 


U 


r2(yi)dy't 


yi-yi> 


dyidy2.  (6.54) 


Pxjiyi  -^Oexpj  -  — 

ii+y2‘ 

Transforming  back  to  the  original  random  variables  X,  and  X3  and  reversing  the  order  of 
integration, 


P*0)  = 


V,  (*(/ 

P\  j(-*3) 

J?i  Jo 


1 


***TZ)» 


rAy’iWi\dxxdx3  (6.55) 


or  if  comparative  fatigue  life  (denoted  in  this  section  by  X  ^)  is  used  instead  of  comparative  crack 
growth  rate, 

/V 


Pdf)  = 


rt(  Co o  (  rxj+t/x  t 

Px,(x  i)  Px*(**)exp  x*  r1(y[)dy[ 

J  t,  J  H(t,-X3)  l  J*3  J 

The  corresponding  expression  for  the  risk  of  static  failure  by  fatigue  is. 

Cl,  C  oo  (  [Xi+tix*,  \ 

P»U)r,(t)  =  I  Px,(*3)  I  r2(x3+?/xT)Pxt(xT)exp|_JC*  I  r2(y’\Wi\ 

Jf t  J  m~x 3)  l  J  i 

(6.57) 

which  if  the  exponential  term  is  ignored  is  the  same  as  the  approximate  expression  developed 
by  Diamond  and  Payne. 

Further  details  of  the  initial  crack  model  are  discussed,  as  part  of  the  documentation  for  the 
NERF  computer  program,  by  Mallinson  and  Graham.13 


dx\dx  3.  (6.56) 


dx\dx 3 


7.  A  THREE-PARAMETER  MODEL 

A  natural  extension  to  the  model  described  in  Section  6.8  is  to  include  a  random  variable 
representing  strength  variation.  The  resulting  three-parameter  model  has  sufficient  generality  to 
encompass,  as  special  cases,  ail  the  models  discussed  in  this  report  and  is  presented  here  using 
notation  which  permits  a  simple  tabulation  of  the  special  cases.  All  these  models  can  be  deduced 
from  this  tabulation  leaving  the  procedures  outlined  in  Section  3.6  for  the  development  of  new 
models  as  the  need  arises. 


7.1  AsMunptioM  and  Model  Equations 

7.1.1  Basic  Random  Variables 

The  model  is  assumed  to  depend  on  three  random  variables  two  of  which  (Xt,X2)  are  para¬ 
meters  in  the  equation  for  strength  decay  as  a  function  of  time  (e.g.  equation  (4.9))  for  cracked 
structures.  The  third  random  variable,  X3,  is  a  parameter  affecting  the  relationship  between  the 
strength  of  a  given  structure  and  some  known  median  strength  for  the  population  (e.g.  x2 
in  (4.10)). 

The  basic  random  variables  will  be  specified  with  no  further  precision.  This  allows  the 
three-parameter  model  to  be  defined  in  terms  of  its  transformed  random  variables  while  leaving 
the  transformations  unspecified,  facilitating  reduction  to  a  wider  class  of  less  general  models. 


7.1.2  Transformed  Random  Variables 

The  transformed  random  variables  are:  age  (Y,),  comparative  crack  growth  rate  (Y2)  and 
virgin  strength  (Y3).  They  are  assumed  to  bear  the  following  functional  relationships  with  the 
basic  random  variables. 
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Y,  =  Y.CXj.X,,/),  (7.D 

Y2  =  Yj(XltXj),  (7.2) 

Yj  =  Y3(X3).  (7-3) 

For,  given  time,  the  joint  density  function  for  Yt  and  Y2  can  be  written  in  the  form, 

Py,  .YjO'iO's)  =  PYjtoW.O'ilPs);  (7-4) 


for  uncracked  structures  Y,  and  Yj  are  independent.  The  density  functions  i  Py/PjIps)*  PyjO^) 
and  pyj(y3)  can  be  derived  from  the  density  functions  for  Xf,  X2  and  X3  using  equation  (3.62). 


7.1 3  Initial  Age 

For  any  structure,  the  age  when  /  =  0  is  called  the  initial  age  (yg)  where. 


y0  =  .Mx,.*i,0). 

(7.5) 

It  is  assumed  that  the  time  history  equations  are  defined  so  that  for  given  time  y0  can  be  alterna¬ 
tively  expressed  in  terms  of  y3  and  y2,  viz., 

y0  =  yoiyuyiJ) 

and  yt  in  terms  of  y0  and  y2, 

(7.6) 

yi  =  ytiyoJi’1)- 

For  a  given  population,  initial  age  may  be  limited  so  that. 

(7.7) 

yo.l  <  >0  <  -Vo.2- 

7.1.4  Model  Equations  and  Subspace  Boundaries 

For  structures  in  Du 

(7.8) 

rt  =  r,(y}) 

and 

(7.9) 

*m m  fi<  f>- 

For  cracked  structures  in  D2, 

(7.10) 

t2  =  r2{y34>(yi)) 

and 

(7.11) 

<  y3;  0  <  y2  ^  y2y,  Pm  yi  <  Pi  .a 

where 

(7.12) 

Pi.i  -  max{f„yl(y0.iJ’2)} 

(7.13) 

y,j2  =  min{fr^,y1(y0.j.>’2)} 

and  y22  is  defined  by. 

(7.1 4) 

PlO’O.lO'2,2.0  =  7f.M- 

(7.15) 

Fatigue  life  limiting,  or  the  effects  of  perfect  inspections  are  included  in  the  parameter 
lij,  where. 

=  min{ft,^-1(P_i,/y3)} 

and 

(7.16) 

H  =  min {!,Myi,t,tij)}. 

The  function  lt^y2,t,tij)  represents  the  effect  of  an  inspection  at  t  =  tij. 

(7.17) 
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7.2  An  Integral  Expression  for  Ps(t) 


Using  the  model  equations,  the  following  equations  can  be  derived, 
and 


where 


Pn(y,‘)  =  exp 

| -ViCKs)  ~ 

't 

t j 

=  exp 

[-Vito)- 

\  • 

M-.»  <«*  *‘J 

=  fl,(yid'2.yi,0. 

r,  =  max{0,/:>’i(.r1,*J,/)  =  f,}. 


The  expression  for  the  probability  of  survival  is, 


(7.18) 


(7.19) 

(7.20) 


Ps(D  =  [ 

J 


exp{  -rl(yJ)t)pyi(yi)dyiPyl(li)+ 


IX.  2 

T’vjfyj)  Pva(y2)  Pyi(yi\yi)HJy,t)dyidy2dyi.  (7.21) 

J  JO  J  ytht 


7.3  Risk  Rate  Equation 


The  expressions  for  virgin  risk,  the  risk  of  static  fracture  by  fatigue  and  the  risk  of  fatigue 
life  exhaustion,  obtained  by  differentiating  (7.21)  with  respect  to  time  are  given  by 


Pn(‘)r.(t)  = 

r  i(yj)exp{  — r  i(y3)t  }Py}(y3)dy3Pyl(,li), 

(7.22) 

J  ^mii» 

l 

*«>  f  yi.i 

fy  i.s 

Pu(‘)r,(t)  = 

Pvjys)  Pv2(y2)  r1{yi)pyl(yi\yl)H£i,t)dyldyidyi 

(7-23) 

J 

4.1.  J  o 

J  y  i,i 

and 

Pu(t>At)  = 

'R^mm 

PY,{y3) 

x.2 

PY,(y2>~ 
y2.f  s' 

PY,(f*|y2)^(f„y2,yj,tVy2<fyj-f 

• 

^mln  « 

+#W» 

f  °° 

Py.X) 

(X* 

y2.f 

PY1(ff!y2W7r.y2.yj.'M)’2^y3 

(7.24) 

J  Rmjm 

where 

=  <r‘(*.i„/y2) 

(7.25) 

and  y2,(  is  given  by, 

yi(yo,2.y2^.0 

=  ?f.,v 

(7.26) 

7.4  Reduction  to  Simpler  Models 

Equations  (7.21)  to  (7.25)  define  the  probability  of  survival  and  risk  rates  for  any  three- 
parameter  model  satisfying  the  assumptions  in  Section  7.1.  In  particular,  all  the  models  described 
in  this  report  can  be  considered  as  simplifications  of  this  three-parameter  model. 
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Table  7.1  details  the  functional  forms  required  to  effect  the  simplifications  for  models 
grouped  in  the  following  manner.  The  column  headed  “Ford  (including  extension)”  incor¬ 
porates  all  the  models  described  in  Chapter  6  that  use  crack  growth  rate  and  initiation  time  as 
the  basic  random  variables  associated  with  age.  The  term  “initial  crack"  identifies  models  that 
follow  Diamond  and  Payne3  and  use  initial  age  and  comparative  fatigue  life  as  the  age  random 
variables.  The  third  column  in  the  table  details  the  two-parameter  model  described  in  Chapter  4. 
Hooke’s  model  (Chapter  S),  is  also  covered  by  the  entries  in  this  column  if  R0  is  set  equal  to  1 
and  X3  interpreted  as  virgin  strength. 

Note  that  the  initial  crack  model  includes  the  strength  random  variable  whereas  the  extended 
Ford  model  does  not.  This  follows  the  model  definitions  of  the  originating  authors  and  this 
random  variable  can  be  included  or  neglected  by  using  the  relevant  table  entries  from  either 
model.  In  this  sense,  the  first  two  columns  describe  the  same  model,  the  essential  difference 
between  the  table  entries  being  the  definition  of  the  basic  random  variables. 


TABLE  7.1 

Functional  Forms  Required  to  Reduce  the  Three-parameter  Model  to  the  Extended  Ford,  Initial 
Crack  (after  Diamond  and  Payne)  and  Payne  et  al.  Two-parameter  Models 


Model 

Ford 

Initial 

Payne  et  al. 

(including  extension) 

crack 

(two-parameter) 

x, 

Crack  growth 

Comparative  fatigue 

Comparative  fatigue 

rate 

life 

life 

Initiation  time 


MXa;  Y,  <  f, 
f.+Y^-X,); 
Y,  ^  *i 


Y* 

(crack  growth  rate) 


Y3 

(virgin  strength) 


Pv.O'ill'i ) 


p*t(.y») 


J'oO'iJ'j.') 


yi(yo,yi,>) 


-■Hr8) h 


Px&i) 


n.a. 


Px,(t) 


yi-yi* 


Initial  age 


Relative  residual 
strength 


X2+t;  Yj  <  f, 

Xj+Yjf; 

Y, 


1/X, 


AqXj 


Px£y\  -ytO 


Px,0/yt)/A 


PXiO's/^o)/A> 


yt-yt* 


Relative  residual 
strength 


&oXs 


Pxt(‘/yi)i/yi 


n.a. 


Px,(y»l&o)l&o 


PxM) 


yo+ytt 

.vo+jy 

n.a. 

i*+yi(t-tb) 

h+yi(t-tij) 

'-(yi-Vlyi 

t-iyi-Wyt 

Myi 

Table  7.2  details  the  integration  limits  appearing  in  equations  (7.21)  to  (7.25).  The  extended 
Ford  model  has  been  subdivided  into  a  model  with  no  initial  cracks  (the  original  Ford1 1  model), 
and  one  allowing  for  some  of  the  population  of  structures  to  be  initially  cracked.  The  initial 
crack  model  assumes  all  structures  are  initially  cracked.  Bearing  in  mind  the  comments  in  the 
penultimate  paragraph,  the  first  three  columns  of  Table  7.2  can  be  used  with  either  of  the  first 
two  columns  in  Table  7.1  to  generate  the  expressions  relevant  to  the  desired  degree  of  initial 
cracking.  Note  that  the  effect  of  the  strength  random  variable  is  detailed  in  column  3  for  the 
initial  crack  model  only. 

All  the  entries  in  Table  7.2  are  relevant  for  fatigue  life  limiting  which  can  be  removed  by 
setting  lf  —  oo,  and/or  Rmim  —  0  appropriately. 


TABLE  7.2 

Integration  Limits  Required  to  Reduce  the  Three-parameter  Model  to  the  Extended  Ford,  Initial 
Crack  and  Payne  et  al.  Two-parameter  Models 


Model 

Ford 

(including  extension) 

Initial  crack 

Payne  et  al. 

No  initial 
cracks 

Some  initial 
cracks 

yo.t 

0 

0 

0 

y0.2 

h 

l, 

0 

yi.t 

h 

h 

max{flff,+y2t} 

h 

>1.1 

min{fr,f,+y2/} 

f, 

>1.2 

«t,-Wt 

n.a. 

yr 

hit 

0 

0 

n.a. 

8.  CONCLUSIONS 

The  method  developed  in  Chapter  3  for  the  generation  of  reliability  models  has  been 
demonstrated,  by  example,  to  be  successful  in  producing  correct  forms  of  the  reliability  functions 
relevant  to  models  used  previously  at  ARL.  In  particular,  the  systematic  application  of  a  unified 
method  has  led  to  an  appreciation  of  differences  between  the  models  and  in  some  cases,  correc¬ 
tion  of  previously  published  expressions. 

The  two- parameter  model  developed  in  Chapter  4  using  the  assumptions  of  the  Payne 
et  al '. 1-4  model  was  shown  to  yield  expressions  presented  by  Payne  and  Graham4  for  the  risks 
of  static  failure  by  fatigue  and  fatigue  fracture.  The  latter  risk  was  shown  to  have  two  con¬ 
tributions;  that  arising  from  an  imposed  limit  on  fatigue  life  and  a  contribution  arising  from  a 
prescribed  minimum  strength  value.  The  model  developed  in  Chapter  4  also  accounts  for  failures 
of  uncracked  structures  whose  strengths  have  not  been  reduced  below  their  virgin  values.  An 
expression  for  the  resulting  virgin  risk  was  obtained  and  the  expressions  for  the  other  reliability 
functions  account  for  this  effect.  This  treatment  of  virgin  risk  mirrors  that  of  Yang  and  Trapp.14 

Expressions  for  the  probability  density  for  strength  and  the  failure  density  of  strength  were 
obtained  and  the  latter  was  shown  to  reduce  to  an  approximate  expression  derived  by  Diamond 
and  Payne.1 

The  two-parameter  model  described  by  Hooke,  although  using  a  slightly  different  definition 
for  the  basic  random  variables,  was  shown  to  be  basically  the  same  as  the  Payne  et  al.  model 
with  the  exception  of  the  assumptions  regarding  fatigue  life  limiting.  Although  Hooke’s  equations 
for  the  probability  of  survival  and  the  risk  of  static  failure  by  fatigue  do  not  account  for  any 


46 


imposed  limit  on  fatigue  life,  Hooke*  derives  an  expression  for  the  risk  of  fatigue  life  exhaustion 
which  accounts  only  for  fatigue  life  limiting.  (No  minimum  value  of  strength  is  prescribed.)  In 
this  sense,  Hooke’s  equations  are  inconsistent. 

An  examination  of  Hooke’s  proposal*  to  isolate  the  failures  of  structures  that  are 
weakened  by  fatigue  cracking  and  fail  by  a  load  less  than  their  virgin  strengths  concluded  that 
the  approach  based  on  taking  the  algebraic  difference  between  two  probabilities  of  survival,  one 
with  fatigue  cracking  and  one  without,  overestimated  such  failures.  Correct  isolation  can  be 
obtained  by  defining  the  basic  risk  rates  used  in  the  model  appropriately  and  generating  the 
corresponding  reliability  functions  using  the  techniques  described  in  Chapter  3. 

Using  the  basic  random  variables  defined  by  Ford,9*11  the  equivalence  between  the  method 
described  herein  and  his  technique  based  on  solving  a  partial  differential  continuity  equation 
for  the  probability  density  for  crack  length  was  established.  Transformation  of  the  random 
variables  led  to  forms  of  the  reliability  functions  which  could  be  readily  compared  with  the 
Payne  el  al.  functions.  By  introducing  initial  age  as  a  random  variable.  Ford’s  model  was 
extended  to  account  for  structures  that  are  initially  cracked.  This  extended  model  encompasses 
the  initial  crack  model  suggested  by  Diamond  and  Payne3  for  the  analysis  of  the  reliability  of  a 
population  of  structures  all  of  which  commence  life  cracked. 

The  three-parameter  model  presented  in  Chapter  7  encompasses  all  the  models  described  in 
this  report;  a  tabulation  of  the  simplifications  relevant  to  each  model  was  presented.  For  these 
models,  this  tabulation  fulfills  one  of  the  objectives  of  the  investigation  by  providing  a  basis  for 
the  complete  specification  for  the  numerical  evaluation  of  the  reliability  functions.  For  other 
models,  or  for  variations  (such  as  different  inspection  procedures)  of  the  tabulated  models, 
the  method  described  in  Chapter  3  provides  the  fundamental  means  for  fulfilling  this  objective. 
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NOTATION 


a  or  o(r) 
a' 


0| 

Of 

A 

Cl 

Dk 


D ! 

D2 

D> 

M 

dFk(x,t) 

/(z) 


F(xuyut) 


«(z) 


Crack  length,  function  of  time 
da/dt 

Crack  length  immediately  after  initiation 

Crack  length  corresponding  to  imposed  fatigue  life  limit 

Random  variable  representing  crack  length 

Scaling  constant  for  the  ith  random  variable — used  in  transformations 

kth  subspace  of  sample  space  for  the  random  variables  describing  the  variation 
in  the  reliability  model.  Each  subspace  represents  a  distinct  time  zone  in  the 
fatigue  process.  For  the  models  described  herein,  three  subspaces  are  used 
as  defined  below. 

Subspace  containing  uncracked  structures 
Subspace  containing  cracked  structures 

Subspace  containing  structures  which  have  passed  a  failure  criterion 
Crack  growth  rate  function,  equation  (6.1) 

Ps(*'‘)Px(x)dx,  see  equation  (3.3) 

Function  of  z  used  in  the  derivation  of  the  failure  density  of  a  structural 
characteristic 

Homogeneous  function  defining  the  transformation  from  xt  to  yu  see  equation 
(3-71) 

Function  of  z  used  in  the  derivation  of  the  density  function  for  a  structural 
characteristic 


//(x,f)  (or  H(x ,,  x2,  ■ . .,  xv,  f)).  Loss  factor  for  x,  see  equation  (3.19) 

H,(t)  Heaviside  step  function 

Loss  factor  for  y 

j  or  i  (subscript)  ith  random  variable  or  initiation,  depending  on  context 
IiCkj.O  Inner  integration  over  x„  see  equation  (3.31) 

k  Subscript  or  superscript,  identifies  subspace 

K  Number  of  subspaces  in  model 

/,  Average  frequency  of  load  application 

L  Random  variable  representing  applied  load 

L  Sample  value  for  L 

M  Number  of  random  variables  in  model 

n  (superscript)  Denotes  a  definition  relevant  to  a  finite  population  of  structures 
P{A)  Probability  of  A 

PA(a)  Probability  that  A  <  a 


A(«) 

Pa(o) 

P utty) 

PdO 

MO 

Pt.% Vs) 

PdX) 
/*(R|0 
Pd 0 
Pity. ) 
pity) 
pi/tyo 

PA*) 

Py.O'iI.Kj) 

PYi.Yjfj’ltJ’i) 

pMo 

pz(z|F  >  t) 

r(0 

%(0 

'c 

MO 

'.(0 

MO 

ftf) 

'•jW 

'■i<y) 

'*(v) 

'ifjtO 

r’tyt) 

R 

R 

* 

A 

A 

Ah. 

S 


Probability  that  A  >  a,  (=1  — PA(a)) 

Density  function  for  crack  length 

Probability  of  detection  at  the  inspection  at  t  =  tij 

Probability  of  failure,  (=P(F  <  /)) 

Density  function  for  the  time  to  failure 

Joint  density  function  for  time  of  failure  (F)  and  structural  characteristic  (Z) 

Probability  that  an  applied  load  exceeds  strength  R,  (=P( L  >  R)) 

Failure  density  for  strength 

Probability  of  survival,  (=P(F  >  /)) 

Probability  of  survival  immediately  before  an  inspection 

Probability  of  survival  immediately  after  an  inspection 

(or  Pg(y,0)  Probability  of  survival  for  structures  with  X  =  x  (or  Y  =  y),  in  the 
subspace  Dk 

Joint  density  function  for  the  random  variables,  X 
Conditional  density  for  Y4  given  Ya  =  y2 
Joint  density  for  Yt  and  Y2 
Failure  density  for  the  structural  characteristic,  Z 

Conditional  density  for  Z  given  survival  to  t.  Loosely  referred  to  as  the  density 
forZ. 

Risk  rate  at  time  t 
Additional  risk 
Constant  risk  (additional) 

Risk  of  fatigue  life  exhaustion  at  time  t 
Risk  rate  in  the  £th  subspace,  Dt 
Risk  of  static  failure  by  fatigue  at  time  t 
Risk  of  ultimate  failure,  see  equation  (5.23) 

Virgin  risk  at  time  t 

Risk  rate  for  uncracked  structures  with  X  =  x 
Risk  rate  for  uncracked  structures  with  Y  =  y 
Risk  rate  for  cracked  structures  with  X  =  x 
Risk  rate  for  cracked  structures  with  Y  =  y 

Risk  rate  for  cracked  structures  as  a  result  of  loads  less  than  virgin  strength 

Random  variable  representing  strength 

Strength 

Median  strength 

Virgin  strength 

Median  virgin  strength 

Minimum  value  of  strength 

Random  variable  representing  survival 


S(x,tij) 

t 

fo 

t‘ 

U 

h 

u 

u 

ft* 

n 

h 

fh 

u  (superscript) 
w  (superscript) 

X, 


*< 

*1.1 

*1.2 

X 

x 

Y, 

Yt 

.Vo 

*(>i) 

'  (prime) 


Removal  function  representing  the  effect  of  an  inspection  at  time  tij 
Time 

Initial  time 

Dummy  integration  variable  for  time 
Time  of  transition  from  subspace  Dk  to  Dk+l 
Median  time  corresponding  to  inspection  criterion 
Median  initiation  time 
Fatigue  life  limit  (=t2  below) 

Median  fatigue  life  limit 

Age  limit  dependent  on  parameter  x2.  Similarly 

Time  limit  representing  the  effect  of  a  perfect  inspection 

ft-ft 

Time  of  jth  inspection 

Reliability  function  with  virgin  strength  preserved 

Reliability  function  for  weakened  structures  failing  by  loads  less  than  the  virgin 
strengths 

/th  random  variable  in  model.  Meanings  of  the  particular  variables  change 
between  chapters,  or  models,  but  are  consistent  within  a  chapter 

Sample  value  of  ith  random  variable 

Value  of  x,  at  boundary  between  Z>,  and  D2 

Value  of  x,  at  boundary  between  D2  and  D3 

Vector  of  random  variables 

Vector  of  sample  values  of  x, 

ith  transformed  random  variable 

Random  variable  representing  age 

Initial  age 

Value  of  yt  at  boundary  between  Dk  and  Dk+l 
Median  relative  strength  decay  as  a  function  of  age 

(1)  Used  to  denote  new  or  different  form  of  existing  function 

(2)  Denotes  dummy  integration  variable 

(3)  Denotes  differentiation  with  respect  to  time 


S(x,tij) 

t 

{o 

t' 

U 

h 

f, 

U 

u 

if*  j 

n 

u 

•h 

u  (superscript) 
w  (superscript) 


X, 


*1.1 

*1.2 

X 

x 

Y, 

Y, 

yo 

>ijk 

iW.fi) 

'  (prime) 


Removal  function  representing  the  effect  of  an  inspection  at  time  tij 
Time 

Initial  time 

Dummy  integration  variable  for  time 
Time  of  transition  from  subspace  Dt  to  Dk  +  l 
Median  time  corresponding  to  inspection  criterion 
Median  initiation  time 
Fatigue  life  limit  (=/2  below) 

Median  fatigue  life  limit 

Age  limit  dependent  on  parameter  x2.  Similarly  ff  i„,  itJ1. 

Time  limit  representing  the  effect  of  a  perfect  inspection 

Time  of  yth  inspection 

Reliability  function  with  virgin  strength  preserved 

Reliability  function  for  weakened  structures  failing  by  loads  less  than  the  virgin 
strengths 

<th  random  variable  in  model.  Meanings  of  the  particular  variables  change 
between  chapters,  or  models,  but  are  consistent  within  a  chapter 

Sample  value  of  ith  random  variable 

Value  of  x,  at  boundary  between  Dt  and  D2 

Value  of  x,  at  boundary  between  D2  and  D3 

Vector  of  random  variables 

Vector  of  sample  values  of  x, 

ith  transformed  random  variable 

Random  variable  representing  age 

Initial  age 

Value  of  yt  at  boundary  between  Dk  and  Dk+l 
Median  relative  strength  decay  as  a  function  of  age 

(1)  Used  to  denote  new  or  different  form  of  existing  function 

(2)  Denotes  dummy  integration  variable 

(3)  Denotes  differentiation  with  respect  to  time 
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